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ABSTRACT 

We assemble 121 spectroscopically-confirmed halo carbon stars, drawn from the 
literature, exhibiting measurable variability in the Catalina Surveys. We present their 
periods and amplitudes, which are used to estimate distances from period-luminosity 
relationships. The location of the carbon stars - and their velocities when available - 
allow us to trace the streams of the Sagittarius (Sgr) dwarf spheroidal galaxy. These 
are compared to a canonical numerical simulation of the accretion of Sgr. We find that 
the data match this model well for heliocentric distances of 15-50 kpc, except for a 
virtual lack of carbon stars in the trailing arm just north of the Galactic Plane, and 
there is only tentative evidence of the leading arm south of the Plane. The majority of 
the sample can be attributed to the Sgr accretion. We also find groups of carbon stars 
which are not part of Sgr; most of which are associated with known halo substruc¬ 
tures. A few have no obvious attribution and may indicate new substructure. We find 
evidence that there may be a structure behind the Sgr leading stream apocentre, at 
^100 kpc, and a more distant extension to the Pisces Overdensity also at ~100 kpc. 
We study a further 75 carbon stars for which no good period data could be obtained, 
and for which NIR magnitudes and colours are used to estimate distances. These data 
add support for the features found at distances beyond 100 kpc. 

Key words: stars: carbon - stars: variable: general - Galaxy: halo - galaxies: indi¬ 
vidual: Sgr dSph 


1 INTRODUCTION 


In the standard ACDM cosmology (Bullock, Kravtsov, & 


Weinberg 2001 White & Frenk 19911 large galaxies form 


through the merger and accretion of smaller galaxies. We 
are fortunate to be witness to this process with the ongo¬ 
ing accretion of the Sagittarius (Sgr) dwarf galaxy ( Ibata, 
Gilmore, fc Irwin[1994 Majewski et al.|2003 1 onto the Milky 


Way (MW). Its proximity provides a unique opportunity 
to study accretion in detail, through the tidally stripped 
streams that straddle the sky. 

There have been considerable efforts to create simula- 


20101 

Lokas et al. (20101; Purcell et al. (20111 

Ibata et al. 

20131 

Vera-Ciro & Helmi|(|2013|). However, the model most 


used by the community is that of Law & Majewski (2010a I 


[hereafter LMIO], designed to fit much of the observational 
data known to that date. 


* E-mail: avon@ari.uni-heidelberg.de 


A recurring theme in the study 
is their use as a probe of the mass 
dark matter halo, for knowledge of t 
is invaluable in constraining such r 

of the tidal tails of Sgr 
and shape of the MW 
lie orbit of the Sgr dSph 

nodels (e.g. [Ibata et al. 

2001| Martmez-Delgado et al.||20041 Johnston, Law, & Ma- 

jewski||2005 Carlin et al.||2012 

|Deg & Widrow|20131 Ibata 

et al.|2013 Vera-Ciro & Helmi 

2013 

1. But no consistent and 


agreed solution has been achieved. Further observational 
data on the actual orbit, especially at larger Galactocen- 
tric radii, will provide additional constraints on models of 
its accretion, and hence the dark matter profile of the MW. 


Many stellar tracers have been used to study the Sgr 
Streams. An early study used carbon stars ([Ibata et ah 


20011 but significant progress came with Majewski et al. 


(20031 who found that M-giants revealed details of the lead¬ 
ing and trailing tidal streams (arms) - data that informed 
the LMIO model. Since its publication, the LMIO model has 
been tested against observations of stars - their positions, 
proper motions and velocities. It captures many of the main 
features of the streams, and has become a benchmark for 
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such studies. Indeed, many regions of the streams are now 
relatively well characterised up to a distance of ~50 kpc by 
a variety of other tracers including horizontal branch stars 
i Ruhland et al.||2011 Shi et aL]|2C)12 l, RR Lyrae |Vivas fc 


et al. 20031, and whose original classification was determined 


Zinn 


2006 Drake et al. 2013a I, red clump stars dCorrenti 


et al. 


2010 Carrell, Wilhelm, & Chen 20121, upper main- 


sequence and main-sequence turn-off stars (Belokurov et al. 
2006 Koposov et al.|[2012 Slater et al.|2013|. 

Despite its success in reproducing many features of the 
observed data, the LMIO model has some limitations. For 
example, it does not explain the “bifurcation” seen in the 
actual stream (Belokurov et al?]|2006 Koposov et al.||2012 


Slater et al. 20131, prompting suggestions (Koposov et al. 


2012 Newby et al. 20131 that there may actually be two pro¬ 


genitors, with separate streams emanating from each. Other 
deviations from the LMIO model have also been found: Using 


RR Lyrae, Sesar, Juric, & Ivezic (2011b) find no evidence of 


the leading arm at low right ascension - although expected 


from the model. Similarly, Drake et al. (2013aI, see no sug¬ 


gestion of the Sgr trailing arm in RR Lyrae in the sky just 
north of the Galactic Plane. 

Beyond about 50 kpc, the data are more ambiguous. 


For example, Newberg et al. (2003); Drake et al. (2013bI; 


Belokurov et al. (2014) each describe a halo feature, located 


at a distance of ~90 kpc (the “Gemini Arm”), which has 
been recently attributed to the trailing arm of the Sagittar¬ 
ius stream. If the Gemini Arm were part of the trailing arm 
of the Sgr stream, it will have a great impact on the possible 
orbit of the progenitor. These data are far from the model, 
as seen on the sky and in distance, indicating that the LMIO 
model requires serious revision. 

Glearly there is much work still to be done to under¬ 
stand both the nature of the progenitor of Sgr, and its fateful 
journey around the Galaxy during accretion. 

1.1 Carbon stars 

In this paper we return to the study of carbon stars as 
a tracer population, the first to reveal the extent of the 
streams ( |Ibata et al. 1120011 . Specifically, we focus of those 
carbon stars that exhibit clear variability. There are a num¬ 
ber of reasons for doing so. Firstly, there are newly published 
lists of halo carbon stars, enlarging the sample. Secondly, the 
ability of the period of any variability to give a good esti¬ 
mate of distance is now better understood thanks to studies 
of the Magellanic Clouds. And, finally there is now an ex¬ 
cellent set of ‘almost” all-sky, variable star, data available 
in the Catalina Surveys ( Drake et al.||2013a l. Carbon stars 
possess properties that make them good probes. They have 
photospheric C/0 ratios that exceed one, and so show dis¬ 
tinct features of C 2 and CN in their spectra, making their 
identification simple even in relatively low-resolution spec¬ 
tra. The strong band heads also aid velocity determination 
( Aaronson| 1983 1, making them good kinematic probes. For 
example, they have been used to study the kinematics of 


the outer regions of the Galactic disk (Demers & Battinelli 


by the presence of certain spectroscopic features. N-type 
stars have substantial absorption at the blue end of the spec¬ 
trum, while R-types stars still show the blue end of the spec¬ 
trum, and CH stars have strong CH absorption. The J stars 
exhibit very strong isotopic bands of C 2 and CN. 

There are two types of explanation as to how carbon 
finds itself on the surface of the star: intrinsic and extrinsic. 
The intrinsic carbon stars are formed by internal processes, 
in which the carbon products of helium burning mix into 
the convective envelope and are brought to the surface. One 
major process is the third dredge-up (TDU), which occurs 
on the asymptotic giant branch (AGB), creating the classic 
N-type carbon stars, which typically pulsate and are losing 
mass. Recent work ( Zamora et al.|2009[ Knapp, Pourbaix, fc 
Jorissen]|2001 1 suggests that the cooler, late-type R carbon 


stars are the same as the N-type, although possibly sitting 
nearer to the bottom of the AGB phase. The warmer, early- 
type R stars are, by contrast, near the red clump, and may 
be formed through an anomalous He-flash. 

Extrinsic carbon stars are formed by the accretion of 
carbon-rich material from a donor star within a binary sys¬ 
tem. These are normally identified (in stars at high Galactic 
latitudes) by high radial velocities - GH stars, for example 
- ( McClure|1984 ). Similar are the nearby so-called ’’carbon 
dwarfs” (dGs). The first was discovered by Dahn et al. (19771 
and was subsequently found to be a binary system (Dear- 
born et al.|1986 1, in which mass transfer has taken place on 


to a dwarf star. These dGs can, however, be identified by 
their high proper motion. 

The AGB carbon stars are the main focus of this paper. 
They are luminous and so can be detected to large distances 
( |Grebel|2007| ), and our primary interest is in the outer halo 
of the MW. Due to a minimum stellar mass required for 


the TDU to occur (e.g. Iben & Renzini 1983 Straniero et 
|al.||2003| ), such carbon stars are also seen as representatives 
of young to intermediate-age populations - and hence of 
the accretion of younger stars than is typical of the halo. 
It should be noted, however, there is still some uncertainty 
about the minimum mass with new models allowing carbon 
stars at lower mass, at low metallicity ( Karakas|2010 |. 

Unsurprisingly, given the variety of carbon star types, 
early work using carbon stars to probe the MW halo con¬ 
stantly encountered the problem that they possess a wide 


variety of intrinsic luminosities in a range of passbands (Alk- 
snis et al.|1998 Abia|2008 1. Using AGB (N and late-R type) 


stars alone might be expected to reduce the problem. Thus, 


Totten & Irwin (19981 assumed a single value of Mr = —3.5 


mag. But even N-types have a wide range of absolute visual 
magnitudes as determined from Hipparcos distances (|Alks-| 


nis et al. 1998). Mauron et al. (2004) use a standard Mk, 


calibrated against LMG carbon stars, but with the addi¬ 
tional refinement of including (J — Ka) colour. This approach 
gives better results, as NIR data take the mass-loss that is 
characteristic of N-type carbon stars into account. 

The uncertainties of distance determination using car- 


2007), the Galactic halo (Hartwick & Cowley 1985 Mould 

bon stars can be considerably reduced if the carbon star is 
a long-period variable (LPV) - as most AGB stars are - as 

et al.|1985||Bothun et al.[|1991||Green et al.|1994p, the LMG 

(Graff et al. 2000 Alves & Nelson 2000 Olsen & Massey we can exploit period-luminosity relationships (e.g. Ita et 

2007), and the SMC (Kunkel, Demers, & Irwin|2000l. 

al.|2004bl|Whitelock, Feast, & van Leeuwen|20081|Soszynski 

The population of carbon stars is composed of a com- 

et al.||2009| |Riebel et al. 20101. But carbon stars are rare 


plex ”zoo” of many types and with a variety of origins (Abia - and those with clear and measurable variability more so. 






































































































































































Tracing MW features with carbon-rich LPVs 3 


The challenge is to garner a statistically signihcant sample 
of such variable carbon stars. This paper tries to do so by 
drawing together all carbon stars in the halo that are to be 
found in the literature. 

Being giants, and often suffering signihcant mass-loss, 
the stars of interest are typically very red. For this reason, 
the largest sample we use is drawn from the large list of 
2MASS-selected carbon stars by Mauron and his collabora¬ 
tors, having colours most consistent with being AGB carbon 
stars. We also draw on other samples, but many of these have 
used spectra at the blue end to identify carbon stars, and so 
preferentially select the bluer early-type R stars, CH stars 
and dwarf carbon stars. Using long-period variables also al¬ 
lows us to be more conhdent that we are dealing with AGB 
stars, and we adopt a minimum period of ~80 days for our 
selection. 

With AGB carbon stars, we should be able to trace 
the Sgr dwarf tidal features to large distances. Carbon stars 
that do not lie on the streams may indicate the presence 
of previously undetected wraps of the Sagittarius stream, 
helping constrain its orbit. Alternatively they could point 
to new, or conhrm previously proposed, substructure. 

This paper proceeds as follows: the data we use are 
described in The methods used to derive distances are 
given in il The results are presented in ^ In we study 
those carbon stars for which we do not have period data. 
The results are discussed in and we summarise in © 


2 DATA 


Our sample of carbon stars are drawn from a variety 
of sources in the literature and are all spectroscopically- 
confirmed. They have very different means of detection, se¬ 
lection criteria and spatial coverage. The Galactic disk, be¬ 
ing young, contains many carbon stars, but in this paper we 
are concerned with potential halo structures, hence we focus 
on those at high Galactic latitudes. 

The final selection of carbon stars is divided into two 
catalogues. There are those stars for which we have measur¬ 
able periods from the Catalina Surveys (hereafter called the 
“first sample”). Note that these stars are listed in Table 
and are identified throughout this paper by the prefix “HG” 
followed by a running number. A “second sample” contains 
those carbon stars for which we do not have good variability 
data. These stars are listed in Table and have the prefix 
“HG2-”. 


For our study, we are interested in the AGB stars, the 
N-type and late-R types. Our samples below rely to a large 
extent on the classifications given by the authors of the var¬ 
ious papers. In some cases (such as Totten & Irwin (1998l) 
each star is specifically classified from its spectrum. In oth¬ 
ers (such as the many Mauron group papers), a general com¬ 
ment is made that properties of the sample are consistent 
with them all being AGB stars. And in a few cases (in partic¬ 
ular those carbon stars found during searches for other types 
of stars), there is no indication of which sub-type they belong 
to. However, we use the colours of the stars to mostly assist 
in attributing them to AGB stars. In the few cases where 
the stars lie outside this colour selection, we use available 
spectra or other information to decide whether to include 
them in our sample. 


For our purposes we require that the stars sit on the 
luminosity-period relationships. Note that we accept “can¬ 
didate” LPV carbon stars, and so may have some CH stars 
in our sample, for example, which may mimic low-amplitude 
AGB stars. However, as we are primarily concerned with 
tracing distant halo structures, we are willing to accept a 
few contaminants. One aim of this study is to locate good 
candidate carbon stars that may trace previously unknown 
halo substructure. These stars can be used as targets for 
deep imaging to reveal any features in the fainter, but more 
numerous, tracer populations. Given that such features may 
be found over the whole sky, having such targets is a very 
cost-effective means of seeking new substructure. 


2.1 Mauron et al. sample 


Our primary source, providing just over half of all our vari¬ 
able carbon stars, comes from the lists published by the 


Mauron group ( 

Mauron et al. 2004| Mauron, Kendall, & 

Gigoyan|20051 jM 

auron, Gigoyan, & Kendall|2007b| (Mauron 

2008 

Mauron et al. 

2014 

1. These were selected from candi- 

dates based on 2M7 

ISS (Skrutskie et al. 20061 photometry 


and then spectroscopically confirmed. Those from |Mauron| 


et al. (20041 and Mauron, Kendall, & Gigoyan (2005 I have 


radial velocities determined, while the later papers do not 


give velocities. For Mauron et al. (2004 1 they had access to 


the second incremental data release of 2MASS, only about 
half the sky; their subsequent work, however, exploited the 
all-sky 2MASS data. 

To identify carbon AGB star candidates they employed 
the 2MASS colours of known carbon stars, which lie in a 
known linear feature in the {H — K3)/{J — H) colour plot, 
selecting stars ±0.15 mag from the median line. Initially 
they limit this sample with the cuts on colour (H — Kg) > 0.4 
and {J—H) > 0.95, and a faint magnitude limit {Ks < 13.0). 
They exclude a region either side of the Galactic Plane, only 
going down to |b| > 25°. For later papers in the series they 
enlarged the selection criteria to Ka < 13.5, |b| > 20° and 
(J — Ka) > 1.3 mag. One star in the Mauron sample (#80 
in their numbering scheme) is actually an S-type, and is 
removed from our list. We also exclude their M56, M57 and 
M58, which they note are known “warm carbon stars” (all 
with a (J — Ka) < 1), and not AGB stars, and which they 
observed to provide comparison spectra across a range of 
carbon star types. 

The published errors o n heliocentric velociti es, where 
given, are ±10-12 kms“^ ( Mauron et al.| 2004 k Certain 
stars have entries that require comment: [Mauron, Kendall,] 
I&: Gigoyan| ( |2005[ ) report a typographic error for star HG2- 
51 (M15 in the Mauron numbering system, which uses a 
preceeding ’M’) in their earlier paper ( Mauron et al.|2004 |. 
Its heliocentric radial velocity is -158 and not ±158 kms“^. 
HG75 (M43) has a greater heliocentric velocity (±54 kms“^) 
in Mauron, Kendall, & Gigoyan (20051 co mpared to the 
more recent observation (±15 kms“°) 


of 


Deason et al. 


(20121. Such differences are not totally unexpected, how¬ 
ever. Signihcant differences 30 kms“^ have been found be¬ 
tween the velocities obtained by absorption and emission 


lines and between optical and radio lines (Menzies, Feast, 
fc Whitelock|2006 1. The reasons are not clear, but they rec¬ 


ommend using absorption lines rather than Balmer emission 
lines for velocity determination. In our table we use the Mau- 



























































4 A. P. Huxor and E. K. Grebel 


ron value of HG75 for consistency with the rest of the table 
sample, as the velocity difference does not significantly af¬ 
fect our conclusions. For HG89 and HG92 we obtain radial 
velocities from Lagadec et al. (20101. 


Mauron et al. 120041 believe the majority of the sample 


in that paper to be AGB stars on account of: their colours, 
the presence in half of Hq emission, and the fact that two- 
thirds have spectra (in which the wavelength range permits) 
that rises to the red side. In their later papers they only note 
that the spectra are typical of N-type stars. 

We also include the carbon stars from the First Byu- 
rakan Survey (FBS) for which the Mauron group under¬ 
took additional spectroscopy ([Mauron, Gigoyan, fc Kendall| 


2007b I. For a few stars, the new spectra clearly show them 


to be N-type. 


2.2 Other Sources 

The Mauron sample provides the largest number of carbon 
stars, but a review of the literature reveals the following ad¬ 
ditional samples. Many of these are also duplicated between 
themselves, or with the large Mauron sample. In these cases, 
the reference for the star given in Table either indicates 
the first report of the star, or the reference that provided us 
with the most data. 


Plane. Cruz et al. (20071 and Reid et al. (20081 essentially 
use the same selections as Cruz et al. (20031. 


Kirkpatrick et al. (20081 selected candidates with the 


following constraints: 


|6| > 25° 

- -50° < Dec < -f5°l 

8<K,< 13.6 (4) 

{J-K,) > 1.3 
and {R — Ks) > 6.5 


Where R magnitudes are obtained from USNO- 
That is, they focus on the southern sky. Later in their pro¬ 
gramme, they reduced the magnitude and colour cuts to: 
9.5 < Ks < 13.6 and {R — Ks) > 5.5 mags respectively. 

A similar, but much smaller, 2MASS-based sample of 
carbon stars was published by Liebert et al. (20001, pos¬ 
sessing “very red” NIR colours, and high Galactic latitudes 
(they too were interested in objects in the halo, not the 
disk). These stars are heavily dust enshrouded. Of these, 
there is only a velocity for HG44 (IRAS 08546-1-1732) from 


Groenewegen, Oudmaijer, & Ludwig (19971. They present 


spectra, and note that all are consistent with being N-type. 


2.2.1 Additional 2MASS-selected stars 


A number of carbon stars have also been found serendipi- 
tously during searches for late-type dwarf stars - based on 


their 2MASS colours - by Cruz et al. (2003 20071; Kirk 
jpatrick et al.j ( |2008| ); [Reid et al.j ( |2008f , and for whom car 
bon stars are contaminants. Fortunately, they published the 
status of these. As contaminants, the authors only state that 
they are carbon stars, and give no indication of the sub-class, 
so they may include CH, dC and early-R type carbon stars 
in addition to the late-R and N-types. 

[Cruz et al.j (|2003P use the following selections: 


| 6 | > 10 ° 
(J -Ks)>l 
{J-H)>1-{J- Ks) 
and J ^ 1.5 x (J — Ks) -I- 10.5 


( 1 ) 


With additional cuts in the NIR colour-colour diagram 
of: 


(J - 77) s; 0.8 
for 0.30 ^{H - Ks) ^ 0.35 

and 


( 2 ) 


(J - 77) ^ 1.75 X (77 - Ks) + 01875 
and (J-77) > 1.75 X (77 -ATs)- 0.4750 (3) 

for 0.35 s: (77 - Ks) ^ 1.20 

They obtained spectroscopy from facilities in both the 
northern (Kitt Peak Observatory) and southern (Cerro 
Tololo Inter-American Observatory) hemispheres. Hence, 
their sample covers much of the sky except for the Galactic 


2.2.2 Totten and Irwin 


Totten & Irwin (19981 list a sample of carbon stars that were 


serendipitously discovered in an APM (Automatic Plate 
Measuring) survey for high redshift quasars (which also ap¬ 
pear as very red, faint stellar-like objects). The APM sur¬ 
vey covered the whole of the northern sky away from the 
Galactic Plane(|6| >30°), and much of the southern sky (see 
Figure 1 of jlbata et al.|[2C)01| . Candidates were selected on 
optical colour [Bj — R < 2.5) and magnitude (11 < R < 17) 
and the equivalent POSSl O and E colours and magni¬ 
tudes. They undertook follow-up spectroscopy and give clas¬ 
sifications for their sample, based on these spectra. We ex¬ 
clude those stars which they classify as CH-type, but include 
HG21 (their 0351-1127) which they classify as borderline 
CH/N-type. 

The radial velocities also come from ITotten fc Irwini 


(19981, except for HG60, where a more recent value is given 


in Lagadec et al. (20101. 


2.2.3 SDSS 


A number of authors have searched through the SDSS 
database and generated catalogues of carbon stars, identi¬ 
fied from the SDSS’s own spectroscopy (Margon et al.|2002 


Downes et al.|20M Green|2013 1. Most of these carbon stars 
are rather blue in colour, due to the SDSS passbands, but a 
handful of AGB stars have been picked up. The SDSS cat¬ 
alogues also provide velocities, but - as Green (20131 notes 


- these are derived from fitting to templates that do not 
include carbon stars. Hence we do not use these values. For 
example, in one case (HG45) [Peason et al. (20121 report a 


^ http://tdc- WWW. harvard.edu/catalogs/ua2.html 










































































































Tracing MW features with carbon-rich LPVs 5 


heliocentric velocity of 200 ±10 kms , which is very differ- 
ent from the SDSS catalogue value of -313 kms“^. We do, 


however, adopt the Deason et al. (2012l value. 


We have spectra for all these, and include those that 
clearly have N-type spectra, and one carbon star that is a 


candidate late-R type (see Section 4.2.11. Many of these are 
among the most distant carbon stars and so we discuss these 
later in Section [|L2] 


2 . 2.4 Others 

The few stars that do not appear in these surveys (and 


mostly listed in Alksnis et al. ( 

20011 ), include: HG2-35 

(CGCS 6480) from Moody et al. 

(19971, found serendipi- 


tously during a search for emission line galaxies; and HG90 
(CGCS 3801) from Meusinger fc Brunzendor^ (2001l found 
during a search for AGN in the fields around two Galactic 
globular clusters. JLee fc Chen| ( |2009[ ) provide one, HG28. 
Christlieb et al. ( 2001[ ) created a catalog of 403 faint high 
latitude carbon stars from a region of 6,400 deg^, from the 
Hamburg/ESO (HE) objective prism survey. This covers 
Dec< +2.5° and |6| > 30°. The vast majority of these are 
not variable, and of those few that are, most are found in 
other lists described above. One exception is HG121 (HE 
2319-1534). 



RA (degs) 


Figure 1. Comparison of |Drake et al.||2014[ | sub-sample of LPV 
stars with our comparable sub-sample drawn from the “first sam¬ 
ple” . Stars known to be oxygen-rich (M stars) are also marked. 


2.2.5 Completeness of sample 

It is difficult to be certain about completeness for a sample 
that is assembled from a set of heterogeneous lists. But one 
limitation we have arises from our use of Gatalina Surveys 
data (see j ]2.5[ ) which excludes the regions around the Galac¬ 
tic Plane, Magellanic Clouds, and the Equatorial Poles. 


Mauron et al. (20041 discuss their sample selection in 


some detail, and we can extrapolate their Hndings to esti¬ 
mate how well we recover the AGB carbon star population 
within the Catalina footprint. Using the colour cuts that we 
described above, they have an initial list of ~1200 objects 
in the 2MASS second incremental data release, which covers 
19,600 square degrees of sky. However many of these are, for 
example, distant galaxies. After cleaning this initial list, they 
hnally have ~200 “good” candidates. The Catalina Surveys 
cover 33,000 square degrees, so we would expect ~340 com¬ 
parable carbon star candidates in this region. [Mauron et al.| 
(20041 also report that some 30% of the 97 stars for which 


they actually obtained spectra were carbon-rich stars, the 
remainder being oxygen-rich. Thus naively one might ex¬ 
pect that just over 100 stars will be found in the Catalina 
Surveys. We actually find 121 in our LPV sample, with a few 
dozen other stars that although carbon stars, do not have 
good light curves. Thus we will have found all carbon stars 
that match the Mauron et al. (20041 “good” candidates. Of 


course, there may be many AGB carbon stars that did not 
match their “good” criteria, which are not given, and where 
the criteria may have been too stringent. 

Another measure of the coverage of our sample is to 
compare it to the sample of the LPVs from [Drake et al.j 
(20141. They use Data Release 1 of the Catalina Surveys, 
which only reached down to a declination of -20°, and so 
does not cover the region south of this (nor the region of 
the northern equatorial pole). We create a test sub-sample 


from our LPV (Erst) sample by applying these cuts for sky 
location. We also apply appropriate cuts to the Drake et 
al. sample with [H — Kb) > 0.4, and |b| >20, to make it 
comparable to ours. This leaves us with two sub-samples 
(Fig.g. Of the 87 Drake et al. stars (Hlled circles), which 
include both carbon-rich and oxygen-rich stars, our sample 
of carbon stars (diamonds) matches to 37, some 42%. We 
can compare this match to Mauron et al. (20041 where some 


30% were shown to be carbon stars. So we believe that we 
are recovering most of the carbon stars in the Drake et al. 
sub-sample, the remainder being oxygen-rich. (Those Drake 
et al. stars that are known to be M stars (oxygen-rich), from 
the SIMBAD or Vizier databases, are marked out with open 
squares in Fig. [^) 

However, what is notable is that we have many LPVs 
in our sample that are not found in the Drake et al. sample. 
We believe a number of reasons may explain our greater 
success rate. Firstly, Drake et al. fail to detect many low 
amplitude stars. They compared their sample to the LIN¬ 


EAR (which we describe below in (2.5 I variability survey 
and found that they missed those stars with an amplitude < 
0.1 mag. Secondly, they appear to have only used the folded 
light curve when visually selecting their variable stars. Many 
of our stars have odd folded light-curves, but are more con¬ 
vincing in the full light-curve (e.g HG24, HG92 and HG104). 
Thirdly, the Drake et al. DRl data covered a shorter dura¬ 
tion (to June 2011: about MJD 55,700 compared to DR2 
which extends to ~MJD 56,500). If a light curve is sparse, 
our longer period will reveal periodicity that may be less ob¬ 
vious in DRl (e.g. HG12,HG14, HG91, HG109 and HG119). 
These light curves are shown in the Online Material accom¬ 
panying this paper. However the light-curves and period fits 
for three stars (HGl, HG2 and HG3) are shown in Fig. 
to show typical examples. 

Nonetheless, there are still many stars in our sample and 
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Figure 2. Light curves from the Catalina Sky Survey for the carbon stars HGl, HG2 and HG3. Left: The x-axis represents the Modified 
Julian Date (MJD) from 53400 to 56700 days. The y-axis represents magnitude in the Gatalina system, where each tick indicates an 
integer value, allowing an estimate of the amplitude of the period. Gentre: A folded light curve over two periods, with the period labelled. 
Right: A plot of the power spectrum (scaled to near to the maximum power), from 50 to 1000 days. The chosen period is indicated by a 
red cross, and the value of the power is labelled. 


that have good light curves in Catalina which did not make it 
into the Drake et al. sample. The lead author of that paper 
(AJD) visually inspected a remarkable ~112,000 stars, so 
it is not unexpected that some true LPV stars might slip 
through the process. So, although we recover the expected 
proportion of their sample, there remains some uncertainty 
as to the completeness of that sample, and by extension to 
ours. 

The fact that our input lists of carbon stars are so het¬ 
erogeneous can partially help here. They combine data found 
in both the optical and NIR, from both the northern and 
southern hemispheres etc. Thus selection biases inherent in 
one will be compensated by others. In summary, it is ex¬ 
tremely difficult to have a true measure of completeness for 
our sample. However, given the arguments above we believe 
that we are finding the majority of the variable AGB carbon 
stars in the region of sky covered by the Catalina Surveys. 


2.3 Carbon stars in MW satellites 

We also investigate the variability properties, if any, of 
known carbon stars in satellite galaxies of the Milky Way. 
These galaxies have well-established distances, and can be 
used to check our procedures. The stars used are listed in 
Table H 


2.3.1 Carbon stars in the Sgr dwarf spheroidal (dSph) 
galaxy 


There are 26 previously published spectroscopically- 
confirmed carbon stars in the Sgr main body ( |Whitelock| 
et al.|[l9^ I. We find that three have detectable periods in 


the Catalina Surveys data; those at the highest right ascen¬ 
sion. The others are not in the Catalina Surveys as they 
lie too close to the Galactic plane. We also have two of the 
variable stars listed in Battinelli & Demers (20131, giving a 
total of five known Sgr carbon stars. 

None of the carbon stars from the Sgr main body given 


McDonald et al. (20121 are in the Catalina Surveys, as 


they also are too close to the Galactic plane, and outside its 
footprint. 


2.3.2 Carbon stars in the Fornax dSph galaxy 


We use the |Whitelock et al.| ( |200^ list of carbon stars in the 
Fornax dwarf. We can identify long-period variable (LPV) 
stars more effectively than Whitelock et al. (20091, as the 


Catalina data are in the optical whilst they used NIR pho¬ 
tometry, and the amplitude of variability is known to be 


greater in the optical than in the NIR (Smith et al. 2002 
[Nowotny et al.|[2011[ ). Also they have 15 epochs over three 
years, whilst we have many more epochs over a seven year 
period. Thus we get periods for stars that they do not. Note 
that one star (star #105 in their table 2) is listed as having a 
period of 340 days, but of 270 days in the caption to their Fig 
9. We obtain a period of 280 days, so we assume the shorter 
value was intended. One of the |Whitelock et al.| ( |200'^ stars 
is also in the Mauron sample (M30 in their numbering sys¬ 
tem - our S19). We also include any additional new carbon 
stars from the Fornax dSph, into our input list, given in 


Mauron et al. (20141. 
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2.3.3 Carbon stars in other dSph satellites of the MW 


Azzopardi, Lequeiix, & Westerlund (19851 provide lists of 


carbon stars for the Sculptor, Carina and Leo II dSph galax¬ 
ies. Sculptor carbon stars are also given in [Menzies et al.| 
120111. Certain Mauron stars are already known to be within 


MW dSph galaxies: their star M4 (our S25) is in the Sex¬ 
tans dwarf galaxy, and is discussed more detail below. Their 


M29 (our SI) is in Sculptor, and discussed in Menzies et al. 


(20111 and Sloan et al. (2012|. As with the Fornax dSph 


star from the Mauron sample, we consider them as part of 
the satellite carbon star sample. 


2.4 Photometric properties of our sample 


We can characterise the ensemble properties of the full set of 
our sample, and compare their properties to our examples of 
MW satellite carbon stars. The stars from the MW satellites 
are not included in the “first” sample, but are used only for 
calibration and comparison. 

Figure [^illustrates the consequences of the different se¬ 
lection techniques made for each sample on their location 
in the NIR colour-colour plots. The upper right of the plot 
is populated by 2MASS-selected stars, while those derived 


from optical imaging - such as the Totten & Irwin (19981 


and the SDSS (Green 20131 samples - are found preferen¬ 


tially toward the blue end. Note that many of the previously 
published carbon stars in the MW satellites are found in the 
blue region, due to their discovery on photographic plates. 

Most of the stars are within the area expected for N- 
type stars. A handful, at the blue end of the distribution, 
overlap in the region of {H — Ks)/{J — H) colour where the 
CH-type carbon stars are also expected (e.g. see Figure j^. 
Inspection of the spectra for the SDSS sample found one 
star (HG77) with a G-band feature, usually considered a 
characteristic of a GH-type star. But as we discuss in §4.2. 1| 
below, it is probably a late-type R carbon star. HG107 is 


an outlier, classified as an N-type by Gigoyan et al. (20011 
(so we include it in our sample) but only with objective 
prism spectra, so this star should be treated with caution, 
especially given its location in the NIR colour-colour plot. 

Notably three stars from the satellites also sit in the 
dG region, one each from Fornax (S6), Sculptor (S2) and 
Garina (S24). The Fornax and Sculptor stars seem to have 
reasonable light-curves, the Garina star less so. The Fornax 
and Sculptor stars also have V-band amplitudes of almost 
one mag, while the Garina star is only 0.38 mag. [Zamoraj 
|et al.| ( |2009| ) find a couple of late-R type stars, which are 
essentially the same as N-type, in the region of the Totten 
NIR colour plot, that is near the dG region, and off the N- 
type, near to these MW satellite stars, suggesting that they 
too may be late-R types. 


2.5 Variability data 

Our star variability data come primarily from Data Release 
2 of the Gatalina Surveys (e.g. Drake et al. 2013aI. The 
main focus of these surveys is the discovery and tracking of 
near-Earth asteroids, but their multi-epoch imaging of large 
regions of the sky provide excellent data for the study of 
the time-varying sky, including variable stars. In total these 
cover ~33,000 square degrees of sky within the approximate 



Figure 3. Colour-colour diagram of (J — H) against {H — Kf], 
uncorrected for reddening. The points are the variable carbon 
stars in Table^ The different symbols represent either the source 
of the carbon star, or in the case of the MW satellite carbon stars, 
the host galaxy. The dotted region illustrates the approximate 
region selected by the Mauron group. They defined their cuts 
from known (but unnamed) carbon stars, so we use their stars, 
resulting from those selections to recreate the selection region. 
The dashed region is the colour selection used by |Cruz et al.| 
l |2003| l. A reddening vector for Av'=5 is also shown, using the 
nor mal|Cardelli, Clayton, &: Mathis| | |1989[ l interstellar extinction 
law ([Binney fc Merrifield|199^ ' 


region —75° < Dec < 65° and more than 10-15° from the 
Galactic plane. The main bodies of the LMG and SMG are 
also excluded, due to crowding. The time coverage is greater 
than seven years. The whole survey has three components: 
the Gatalina Sky Survey (GSS), the Southern Sky Survey 
(SSS), and the Mount Lemmon Survey (MLS), which cover 
the northern sky, southern sky and the ecliptic (±10°) re¬ 
spectively. In many cases, stars have data from more than 
one of these surveys, and in a few cases all three, adding to 
the extent and quality of the light curves. 

The Gatalina Surveys data ( Drake et al.|200^ 2013a|b ( 
are provided by three telescopes: 1) The 0.68 metre Gatalina 
Schmidt Telescope, which has a 8 deg^ field of view, a pixel 
scale of 2.5 arcsec and which reaches a limit of V ~ 19 — 20 
mag. 2) The 1.5 metre Mount Lemmon reflector, with a 1.1 
deg^ field of view goes deeper, to V ~ 21.5 mag. And 3) 
the 0.5 metre Uppsala Schmidt at Siding Spring, Australia, 
with a 4.2 deg^ field of view a 1.8 arcsec/pixel image scale 
and a depth of U ~ 19 mag. 

All three telescopes operate with cooled, 4kx4k back- 
illuminated unfiltered GGD cameras and use comparable 
data reduction pipelines. The telescopes observe for 21-24 
nights per lunation (avoiding bright nights). A sequence of 
four exposures of each field are taken, each separated by 10 
minutes. The sequences are repeated, giving a typical 2 to 4 
times, per field, during each lunation. The Gatalina Surveys 
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Figure 4. Carbon star colour-colour diagram of (J — H) against 
(H — Ks), uncorrected for reddening. Symbols are the same as 
Fig|^ The regions show the locations of N-type, CH-type and dC 
carbon stars, derived fro m Figure 3 of |Totten, Irwin, &: Whit elock| 
| |2000[ |. We see that the [Totten, Irwin, &c Whitelock] | |2000[ | stars 
fit within their region defining the N-types, yet the Mauron group 
stars also fall just below. Also shown is the location of {J — Kb) = 
1.2 (red dotted line), which we use to distinguish AGB stars from 
our “second sample” (in ij^) 


photometry is undertaken with SExtractoi]^ and the data 
are made available through their websit^ 

In one case we also use LINEAR data in addition to the 
Catalina Surveys. The LINEAR survey ( |Sesar et al.||20lla l 
is similar to Catalina, in that its primary aim is the search 
for Near Earth Objects (NEOs). It operates on two 1 metre 
telescopes located at the White Sands Missile Range in New 
Mexico. The camera is unfiltered and has a ~2 deg^ field of 
view and 2.25 arcsec/pixel image scale. However, it is not 
as deep as Catalina, reaching down to E < 18, and the 
current public LINEAFij^ database only covers the overlap 
with the main contignous Legacy Survey area of the SDSS 
DR7 footprint. The only example of the use of LINEAR is 
for HG53, for which the Catalina data conld not well fit (see 
Appendix). The LINEAR fit is comparable to our Catalina 
data fit, with a period of 128 days and an amplitude of 0.58 
mags (Fig. [^. 


^ The Catalina data are unfiltered, but a conversion to the 
Johnson-Cousins E-band is provided by V = Vqss + 1-07 X (E — 
7)2 -I- 0.04(cr = 0.063) 

^ http://nesssi.cacr.caltech.edu/DataRelease/ 

^ https://astroweb.lanl.gov/lineardb/ 
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Figure 5. Light curve from the LINEAR survey of HG53. 


3 DISTANCE DETERMINATION 


3.1 Period determination 

We use the light curves for our sample of carbon stars to de¬ 
rive their periods, from which we determine luminosities and 
hence distances. Light curve fitting was undertaken with the 
PeriodOl code ( Lenz fc Breger||2004 1, which has previously 
been successfully used on the light curves of AGB stars (e.g. 
Kamath et al.pOlOl. The light curves, folded light curves. 


and power spectra for the MW satellite carbon stars, and 
out first sample, are shown in the figures in the Appendix 
to this paper. 

Many carbon stars - estimated at approximately one 
third of carbon Miras, and an unknown proportion of other 
carbon-rich variable stars - exhibit major obscuration events 
i Whitelock et al(1|2006 l. Soszyriski et al. (20111 present ex¬ 
ample light curves for both O and C-rich stars. The C-rich 
stars are seen to exhibit much greater changes in their mean 
luminosity, connected with episodes of intensive mass-loss. 
During these events, the magnitude can drop dramatically. 


Feast, Whitelock, & Marang (20031 conclude that these are 


due to the formation of dust clouds along the line of sight; 
an idea supported by the simulations of|Woitke & Niccolinij 
( [2005] ). 

To account for these obscuration events, we visually in¬ 
spect all our light curve fits. In a few cases, only parts of the 
light curve are used for fitting, as obscuration events make a 
solution for the whole duration implausible. One example is 
HG4, which has a sudden drop of ~ 1.5 mags in mean magni¬ 
tudes at a MJD of ~54830 days, with a slow rise thereafter. 
Hence, in this case, we only fit the data up to that date. The 
drop is so sudden that no long secondary fit will fit all the 
data. 

It is also notable that many of our light curves, from 
spectroscopically-confirmed carbon stars, have regular peri¬ 
ods (e.g. HG46). In jSoszynski et ah] ( |20lf| such regularity 
is seen as characteristic of 0-rich stars. Although they pri¬ 
marily use the stellar colours to assess chemistry, |SoszynsH| 
jet al.j ( |2011[ ) also use light-curve morphology to assign stars 
to either O-rich or C-rich classes. The Catalina light curve 
data suggest that one should be cautious in assuming that 
well-behaved light curves indicate an O-rich star. 


Mauron et al. (20141 list their own period determina¬ 


tions from Catalina data (combined with LINEAR data in 
some instances), of which 63 are shared with our own “first 
sample”. There is an excellent agreement between our peri¬ 
ods and theirs, despite the different analysis techniques used 

(Fig.§. 

An exception is the star HG52 (M7) for which we have a 
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Figure 6. Plot of period from |Mauron et al.|||2014[| (filled circles), 
[Battinelli &c Demers| ) [2012| l (diamonds) and |Battinelli &: Demer^ 
( |2013[ l (squares) against ours, where objects are common. The 
dashed line shows where the periods are equal, and is plotted to 
guide the eye. 


very much shorter period, 166 days (compared to their 314.8 
days) on top of a longer secondary period (LSP). Inspection 
of our fit to the light curve, compared to that of|Mauron et| 


al. (20141 indicates that the shorter period is a better fit to 


the data. 

Another outlier is HG27 (M69) where again we have a 
period about half that found by |Mauron et ah (2014|. We 
believe that the light curve, folded light curve and power 
spectrum (see Appendix) support our shorter period. 

There may be artefacts in the period determination for 
both |Mauron et ar] ( |2014| ) and ourselves, due to the cadence 
of the Catalina observations that we both use, so it is also 
valuable to compare our periods to other data sets. A few 
stars are also shared between our first sample and that of 


Battinelli & Demers (2012 20131, and allow us to compare 


the periods determined. These data allow a good cross-check 
as they were taken in NIR filters, compared to the optical 
band of the Catalina Surveys, and used different light-curve 
analysis tools. Moreover, their data were observed within the 
dates (MJD from 55,250 to 55,950) contained within those 
of the Catalina Surveys data, ensuring a fair comparison. 

Again we see a good match except for a couple of out¬ 
liers, HG70 (Mil) and HG98 (M18). Inspection of the |Bat- 


tinelli & Demers (2012) light curve shows that their data for 


HG70 (Mil) are the least reliable, with few data points, so 
such a bad match to our period is not unexpected. Similarly, 
HG98 (M18), has fewer data points than most of their sam¬ 
ple, and also reaches down to the faintest magnitudes. Our 
light curves show our fit is very good for HG98, but more 
uncertain for HG70. 


3.1.1 Period-luminosity relations in LPVs 

We determine the distances to our “first sample” stars from 
their periods. The primary period of a long-period vari¬ 
able (LPV) star is a good luminosity, and hence distance, 
indicator, although it is known that a number of period- 
luminosity relations exist (e.g. Ita et al.|2004a Fraser, Haw- 
ley, fc Cook|2008 1. Each relation follows a sequence in a plot 
of log(period) against magnitude (see top panel Figure [7|. 


Using the notation of Soszyhski et al. (20051 and Ita 


et al. (2004al there are three sequences of relevance, where 


the period is greater than about 90 days: C', C and D. The 
different sequences in the PL relationship are due to a variety 
of reasons. In particular the C and C' sequences are thought 
to represent the fundamental and hrst overtone of pulsation, 
respectively. The origins of the D sequence are still debated, 
but may involve binarity ( Soszyhski fc Udalski|[20l4 (. 

Our first task is to allocate our carbon stars to one 
of the sequences in the P-L relationship. We achieve this 
by using a large sample of carbon stars with relatively well- 
known properties (and a known distance) from which we can 
derive a likely connection between the observable properties 
of the carbon star (i.e. colour, period, amplitude) and the 
sequence to which it belongs. Overtone variables (C'), for 
instance, are expected to have smaller amplitudes than the 
fundamental (C) mode ( Wood &: Sebo|1996 |. 

Our main comparison set of LMC carbon stars is the 


sample of Soszyhski et al. (20051 from the OGLE survey. 


which they identify from their position in the Period-W/ 
plot, where Wi is a Wesenheit Index (Madore||1982 1, given 
by: 


Wi=I- 1.55 X (U - /) (5) 

This is a reddening-free quantity that is dehned for any 
two given wavebands Al, A2 with magnitudes as 


Wai,a2 — rn\2 — — m\2) 


( 6 ) 


They find that the location of spectroscopically- 
confirmed carbon-rich and oxygen-rich AGB stars separate 
effectively in these plots ( Soszyhski et al.|2005 (. We use their 
sample of LMC carbon stars, cross-matched with 2MASS to 
obtain NIR photometry, to create a series of selections that 
allow us to assign our sample to certain sequences in the 
period-luminosity plot. We remind the reader that the also 
LMC contains much younger populations than expected for 
the Sgr dSph. Hence we cross-check and rehne these cuts by 
ensuring that they attribute C and C' to the MW satellite 
carbon stars such that they give consistent distances. 

Previous applications of the LMC P-L relations have 
mostly used the A-band magnitudes. We, on the other hand, 
follow-up the suggestion of Soszyhski et al. (20071, who 


noted that a Wesenheit index in the NIR allows for the cor¬ 
rection of the reddening due to the circumstellar dust. This 
is because in the 2MASS colours, it is hard to distinguish 


between interstellar and circumstellar extinction (Messineo 


et al. 20051. Thus, we create a Wesenheit Index for J and 


Aa,denoted {Wjk)’. 


WjK = Aa - 0.686 X (J - Aa)o (7) 
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Figure 7. Plots of periods against luminosity for LMC carbon 
stars, from [Soszynski et al^ | |2009[ l data. Top panel: Plot of pe¬ 
riod against Ks magnitude. Centre panel: As above, but using 
the bolometric magnitude, derived from J and Ks, as given by 
[Whitelock et ar] l |2006| l. Bottom panel: As above, but using the 
Wesenheit Wjk magnitude. The blue dotted lines separate, by 
eye, the three sequences: C', C and D. The solid green lines are the 
fits for the sequences C', Wjk = —4.80(/og[per2oci] — 2.0)-t-10.39, 
and C, Wjk = —5.19{log[period] — 2.0) + 12.01 from |Soszyhski| 
|et al.||20^ . 


Using this value we find that the scatter of the LMC 
carbon star sequences, shown in the log{period) — Ks plot, 
tighten up considerably, when we replace Ks with Wjk (see 
bottom panel of Figure]^. This allows for a better discrimi¬ 
nation of stars between the various sequences, assuming that 
the scatter in the sequences are mainly due to circumstellar 
dust. 

We also investigated the use of the relations between 
period and bolometric magnitude (Whitelock et al. 20061. 
We plot the LMC PL relationships using their approach and 
we see that it is not as tight as the Wesenheit index (see 
centre panel of F ig. [7|. We believe that this is because the 
relation given by Whitelock et al. (20061 requires averaged 
values for J and K, which they obtained from multi-epoch 
monitoring for their sample. However, as (J — K) varies 
throughout the phase of pulsation, our single-epoch 2MASS 
photometry will likely have the wrong NIR colour to put 
into the relation. Hence we use the Wesenheit Index for this 
study, as it is less sensitive to these changes in colour. 

If the 2MASS observations occurred during the Catalina 
Surveys, we might be able to estimate where in the light 
curve they were made. Unfortunately, the 2MASS data pre¬ 
cede the Catalina Surveys by many years, making any such 
estimate impossible. 


3.2 Determination of P-L sequence for our sample 

We determine the observational features from the LMC data 
that we can use to assign a sequence to our first sample. This 
is done step-wise in the following manner. 

Firstly, it is relatively simple to separate out the LMC 
D sequence stars as they primarily fall below the line (see 
Figure : 


(J — Ks) = 3.5 X log{period) — 7.8 


( 8 ) 


The few apparent D stars sitting above this line are 
almost certainly C sequence, as they have a large amplitude, 
characteristic of the C sequence. One can see a number of 
very reddened stars in the bottom panel of Fig. which still 
fall below the line separating the C and D sequences. Those 
stars above the line log(period) > 2.4, i.e. greater than ~250 
days, are assigned to the C sequence. 

For stars with a period below 250 days, it is found that 
many of those with I-band amplitudes > 0.39 mags are ac¬ 
tually on the C sequence, as are those with (J — Ks) > 2.0 
mag. This value for the amplitude was adjusted to 0.42 mags 
to ensure that all the calibration stars in our MW satellite 
sample (Table have an appropriate distance. This ap¬ 
proach also allows us to take into account that the LMC 
data have amplitudes in the 7-band, while our data are in 
the V-band. We cannot, therefore, simply apply the selec¬ 
tion that is most appropriate for /-band data. Moreover, any 
correction from I to V-band is almost certainly a function of 
colour. Using our MW satellite data circumvents this issue 
to some extent, although we acknowledge that the allocation 
of sequence adds uncertainties to our results. 

For those that remain, the plot of {J — Ks) against 
log(period) reveals that many of stars are on the C sequence 


if they are below the line given by (Fig. 101: 
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Figure 8. Plot of periods against (J — Ks) mag for LMC carbon 
stars, from [Soszynski et al.|p009[ | data. The green points repre¬ 
sent the sequence, black points are stars on the C sequence and 
blue points are D sequence, as derived from the selections shown 
in Figure]^ The region occupied by sequence D stars is primarily 
below the dotted line: (J — K^,) = 3.5 x log(Period) — 7.9, and 
stars above this line and with a log(period) > 2.4 (dashed line) 
are assigned to the C sequence. 
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Figure 9. Plot of periods against the I-band amplitude for LMC 
carbon stars, with data and symbols as in Figure The stars 
with log(period) <2.4 (dashed line) and an amplitude < 0.39 are 
mainly on the sequence. 


(J — Ks) = (0.9 X period/l^Q) -h 0.4 (9) 

This last selection is the most uncertain. If a star fails 
these cuts, it is assigned to the sequence. Inspection of 
Figure shows significant scatter of C and stars on 
either side of this line. Hence we expect that there will be 
some mis-allocation of C and sequence stars by our cuts 
for these short period LPVs. 



Figure 10. Plot of log(period) against J-K mag for stars with 
periods < 250 days, with data and symbols as in Figure]^ Al¬ 
though with some scatter, stars below the black dotted line (given 
hy J — K = 0.9 X period/-\- 0.4) are mostly C sequence. Also 
shown are the first sample stars that lie in this range of periods 
(red asterisks). 



Figure 11. Distribution of C (black hatching) and (green 
hatching) stars, with data and colours as in Figure separated 
into those found above (top panel) and below (bottom panel) the 
line defined by 0.9 x period/IbO + 0.4. The y-axes are both set to 
140 counts to ease comparison. 

In Figure El we plot the number of LMC C sequence 
and C' sequence carbon stars that are found above and be¬ 
low our final selection cut. We see that for a period < 210 
days, the cut selects the C' sequence relatively effectively. 
Above 210 days, the selection of C sequence stars is also 
reasonable. However, there is considerable mixing of C and 
C' stars below 210 days. 

It is notable that our “first sample” contains very few 
stars on the C' sequence (see Table , compared to those 
in the LMC (se e Figure [7|. Th e study of LPVs in NGC 147 
and NGC 185 (Lorenz et al. 2011| ) reveals that these two 
dwarf elliptical galaxies also contain few C' stars compared 
to the C sequence, which they put down to their age. How- 
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ever, crowding probably affected these observations, which 
is more likely to impact on the smaller amplitude C' stars 
in these relatively distant galaxies. Similarly, our Catalina 
data are also biased against low amplitude stars, where the 
photometric errors become proportionally large as compared 
to the amplitude. 

We do not use the period-luminosity relationship based 
on these LMC carbon stars. Instead, we use the carbon stars 
of our MW dSph satellite sample to rederive the relation. 

The period-luminosity relations in the literature are 
typically based on the LMC. But the LMC experienced an 
increase in star formation ~ 3.5 Gyr ago ( Weisz et al.|2013 l, 
whilst the primary population of Sgr is estimated to have an 
age of ~8 Gyr ( [Bellazzini et al.|2006b |. Moreover, the stars 
in the Sgr stream (one of our main objects of investigation) 
are older and more metal-poor than those in the bound rem¬ 
nant of Sgr ( Bellazzini et al.|2006a Chou et al.|2007 |. Thus, 
we might expect that a calibration based on the LMC may 
not be appropriate. By using a range of other MW satellites, 
we provide a sample with a wider range of star-formation 
histories. Moreover, as there are Catalina Survey data for 
32 LPV carbon stars in these other satellites, compared to 
only 2 or 3 which are likely in the outer regions of the LMC, 
we can use a calibration set that is more comparable with 
our “first” sample, having undergone the same procedures 
as our first sample. 

Using our MW satellite stars requires NIR photometry 


for them, which was obtained from 2MASS (Skrutskie et al. 


2006), except for the two Leo II stars. For these, we use 
the NIR photometry of |GulIieuszik et al. (2008) as, at the 
distance of Leo II, the 2MASS photometry is not reliable. 
The photometry is corrected for extinction and reddening by 
dust in the MW using Schlegel, Finkbeiner, & Davis (19981. 
As our sample lies away from the Galactic plane, these cor¬ 
rections are mostly small. 

We then fit the Wesenheit magnitude (Wjk) against 
the periods found for these stars from Catalina (bottom 
panel of Fig. 12), and determine the best linear fit for the 
C sequence stars (solid line). As only one star sits on the C' 
sequence, we simply add 1.53 mags to the C sequence fit. 
For comparison we also plot (dashed lines) the fits for the 
LMC data as found in Figure Note that the PL relations 
derived from the LMC and the other MW satellites are simi¬ 
lar, with a 0.1 mag offset for a period of 300 days, increasing 
to 0.2 mag for a short period of 125 days. These values are 
typically less than the uncertainties driven by the impact of 
variability. 


3.3 Distance determination and uncertainties on 
distances 

From these absolute Wesenheit magnitudes, we obtain the 
distances for all stars in our first (LPV) sample, which are 
given in the Table 

One major source of uncertainty in the absolute magni¬ 
tude, and thus distance, is the variability of the stars. The 
PL relationships use mean K-band magnitudes. As we do 
not have information on the phase in which the 2MASS ob¬ 
servations were made, they may (in the worst case) be at 
the maximum or minimum of the variability. 

To estimate the variability in the Ag-band amplitude 
of our sample we use their (J — K) colour. This is possible 



log(period) [days] 


Figure 12. Period-luminosity plots in both the Ag-band (top) 
and Wjfc (bottom) for the satellites of the MW for the C sequence 
— excluding the star S7, in Fornax, which has the shortest period 
and lies on the C' sequence. Top panel: Best linear fit (solid line) 
with estimated uncertainties in M/Cg. Bottom panel: The dashed 
lines show the same fits to the LMC data as in Figure The 
solid lines show the best linear fit for the stars on the C sequence 
(IFjif = —4.771 X log[period] +2.698) - the one sigma uncertainty 
on the slope of this fit is 0.42, and of 1.01 mag for the intercept. 
The upper line has the same slope as the lower, but is 1.53 mags 
brighter (the difference for the LMC data at log[period] = 2.3: 
see bottom panel of Fig.[^. 


as there is a correlation between the amplitude of A-band 
variability. We use the data from Table 3 of|Whitelock et| 
[ah] ( [200^ - who had multi-epoch observations covering the 
whole period of variability - for all their stars with ampli¬ 
tudes given (excluding those stars for which only limits are 
given). A linear fit (Fig. |13[) to these data give the following 
relation: 
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Figure 13. Plot of mean ( J — iCs)o colour against K-band ampli¬ 
tude for the |Whitelock et al.| l [^06| | sample (asterisks), with the 
best linear fit (solid line). Also shown is the smaller [Battinelli &:| 
[Demers] ( |2013| l sample (squares) with the best fit to these data 
(dotted line). 


Kampiitude = 0.334 X (J — Ks)q — 0.116 


( 10 ) 


We can check this relation using the sample of |Battinelli| 
|fc Demers] | |2013[ ), who obtained J and K light curves for 
many of the Mauron group stars, and so overlap our hrst 
sample. We see that the fits are comparable, although both 
show some considerable scatter (Fig. |13| ). It is difficult to 
estimate the impact of this scatter on our distances as it 
appears irregular; for example, there is a smaller scatter at 
{J—Ka) ~ 3.7 than either blueward or redward of this value. 
The distance errors are then derived by taking half of the 
peak-to-peak amplitude in K magnitudes and propagating 
these through the distance estimates. 

The plot of period against distance (Figure [Id} bottom 
panel) shows that the more distant halo carbon stars in our 
hrst sample have shorter periods. It has been suggested that 


the period of an LPV is an indicator of age (Menzies et al. 


2011). If this is correct, the more distant stars are older 


than those nearer to the MW. One object stands out here, 
HG118, which has a weak light curve, and which we discuss 
later (if as one of our least certain stars. 


4 RESULTS 

The main results of this study can be found in Table[^where 
details of the variable carbon stars are given. This includes 
RA and Dec (J2000), derived NIR values (from 2MASS), 
our distance estimate, and any radial velocity (converted to 
Galactic Standard of Rest) as given in the literature. We also 
give the period and amplitude (in Gatalina V magnitudes) 
found from our best estimate using PeriodOd hts, and the 
adopted sequence (C or C'). Finally, we also give the refer¬ 



Dlstonce (kpc) 


Figure 14. Top: Plot of {J — Kb)o colour against distance for our 
“first” sample of our LPV carbon stars (Tables 2, 3 and 4; black 
symbols). Bottom: Plot of period against distance. The distances 
shown for the historical carbon stars in known dwarf galaxies 
(open symbols) are those derived by our method, and not of the 
known distance to the host, to illustrate the spread of values 
and systematic errors in our method. The distances, and their 
uncertainties, to the host galaxies l |McConnachie|2012[ ) are shown 
by the same symbols with error bars at the bottom of each panel. 


ence in the literature that provided the most valuable data, 
and note the halo feature to which we judge the carbon star 
to most likely be associated. Note that Unc in this table 
represents “unclear”, where it is not clear to which - if any 
- structure in the star may sit (a colon after the feature 
name indicates, in the usual manner, some uncertainly in 
the attribution given). 


4.1 Comparison of the data to the LMIO model 

One of the major goals of this study is to investigate whether 
the distances derived from newly available Gatalina light 
curves will better reveal those carbon stars that belong to 
the Sgr streams, and those that may be from other substruc- 
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ture. In many cases we have been able to use the location, 
distance and (where available) the radial velocity of a carbon 
star to assign it to known substructure. These assignments 
are given in the final columns of Table For many sub¬ 
structures, they are also indicated in Figure except for 
the overlapping set of candidate structures near the Galactic 
Anti-centre. 


Previous authors, such as Mauron et al. (|2004 1 have 
compared their halo carbon stars to the older |Ibata et al.| 
(20011 model, and later papers compare their carbon star 
data to the model of Law, Johnston, & Majewski (20051. 
For the first time we study a large sample of LPV carbon 
stars comparing it to the LMIO model. 


Plots of the location of the carbon stars in comparison 
to the LMIO model, for location of the sky, distance and ra¬ 
dial velocity, are shown in Figure For these plots we use 
a simple selection to identify stars from our lists that may 
be associated with the disk. Using the distance, we derive 
the locations of our stars in the Galactocentric coordinate 
system. We then highlight stars that are within a radius of 
18 kpc from the Galactic Centre and a distance of 5 kpc 
above or below the Galactic Plane. These are indicated as 
“Disk” in the tables, and by asterisk symbols in Figure [T^ 
Our values for the size of the Galactic disc region is a con¬ 


servative one, based on a disc radius of ~15 kpc ( Ruphy et 
|al.||1996] [Minniti et al.|[2011[), but taking into account the 


known presence of young stars out to 20 kpc (Carraro et al. 


2010). Heliocentric radial velocities, as given in the litera¬ 


ture, were also converted to the Galactic Standard of Rest 
(Vgsr)- 


A notable aspect of the top panel of Figure [T^ is how 
many stars of the first sample do not have published veloci¬ 
ties. Yet we can see how well the velocity data can help iden¬ 
tify the stars that are most likely in the Sgr tidal streams. 
This immediately suggests that a programme of observations 
to obtain radial velocities for all the carbon stars without 
such data will be invaluable. Given the luminosity of these 
stars, such an effort will not be too demanding. 


One significant result is that although the majority of 
our input lists were based on all-sky surveys (e.g. 2MASS), 
the majority of the carbon stars in our study appear as 
though they can be attributed to the Sgr streams. This is 
particularly true for those for which velocity (Vgsr) data 
are available, but is also the case for central (RA against 
heliocentric distance) and bottom (RA against Dec) panels. 
Many carbon stars trace the leading arm across much of its 
length, and a group of them at RA ~300° are very close the 
main body of Sgr (which itself lies near the Galactic Plane 
and hence is not in the Catalina Surveys). In many cases we 
can allocate the stars to either the leading or trailing arms 
(green and red points respectively in Figure [l^ . 

The apocentre of the trailing arm in the LMIO model, 
centred at RA ~100°, contains no carbon stars. This is not 
totally unexpected as much of it lies behind the Galactic 
Plane. But, even taking this in consideration, there are very 
few stars in the region of the trailing arm apocentre just 
north of the Galactic plane. This may suggest that the LMIO 
model is not a good representation of the trailing arm. Note 
that all of the (three) carbon stars that appear to sit in the 
arm at the appropriate distance (middle panel of Figure [T^, 
and that are marked with green open circles, are actually in 


the southern Galactic hemisphere, and associated with the 
LMC (see lower panel). 

It is worth commenting on a number of stars. Two stars 
(HGll and HG23) stand out as inconsistent with LMIO, 
sitting either side RA = 50°, and at a distance of ~60 kpc. 
Both these have properties (period, colour and amplitude) 
that place them close to the boundary between the G and 
C' sequence in Figure 10 They also have Vgsr values (- 


86 and -156 km/s), however, that are consistent with them 
sitting in the trailing Sgr arm. Hence, we believe that these 
are actually C sequence stars, not the C' that our selection 
attributes to them. If this is the case, it demonstrates the 
caution that we must exercise in using these selections near 
the boundary. HG107, which was discussed in §2.4| at RA~ 
320° and a distance of 50 kpc is quite likely (given its NIR 
colours) not an AGB star. 

We can also look at the locations of our first sam¬ 
ple stars in Galactocentric cartesian coordinates (Fig. |16|. 
Again it is striking that the majority of our sample is con¬ 
sistent with much of the LMIO model. The distant carbon 
stars that we attribute to A16, the Gemini Arm are at the 
edge of the plots. It is also easier, in this representation, 
to see the apparent lack of Sgr stars - in the X vs. Z plot, 
just above the Galactic Plane and affecting both the lead¬ 
ing and trailing arms. The two stars that we believe may be 
mis-assigned G' can also be seen. 


4.1.1 The Nature of Distant Carbon Stars 


There are six carbon stars with estimated distances that are 
greater than 100 kpc. Some of these are somewhat blue com¬ 
pared to the rest of the first sample, which arouses suspicion 
that they may be CH, early-R or dG stars. 

The spectra for four of these six stars are available from 
the SDSS, and we also plot the spectrum for another SDSS 


star at 68 kpc (Fig 171. Four show classic N-type spectra, 
with little flux below 4500A, and all four of the stars more 
distant than 100 kpc showing Balmer emission. All four ex¬ 
cept HG66, show very good light curves. Although the spec¬ 
trum of HG45 has a low signal-to-noise, it is consistent with 
a carbon star. Its light curve (see Appendix) is also of good 
quality, reinforcing its claim to be an AGB star. 

The distant carbon stars are apparently associated with 
other tracers, and all are likely related to structures which 
are already known. These are labelled on the middle panel 
of Figure [T^ as 1, 2 and 3. Many other stars of the sample 
can be attributed to known substructures. We now discuss 
each of these groups and features in more detail, followed by 
the closer substructures. 


4 . 1.2 A16/Pisces Overdensity (Feature 1) 


Feature 1 comprises two carbon stars at a distance of ~100 
kpc at RA ~ 10-30. These are probably associated with a 
feature known as A16, detected in M-giants by |Sharma et al.| 
(20101, who show it to be a rather large and diffuse structure 
at a distance of 100 kpc and stretching 20° either side of RA 
= 0. A16 encloses the previously reported and better known 


“Pisces Overdensity”, a feature found in RR Lyrae (Watkins 
|et al.|[^09[ |Sesar et al.|2010b| ) in the SDSS stripe 82, lying 
at a distance of ~80 kpc, and by Drake et al. (2013b I in 
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Figure 15. Location of carbon LPVs in relation to the LMIO model (small points) for the last wrap; showing the leading (green) and 
trailing (red) arms. Top panel; RA against the radial velocity (with respect to the Galactic Standard of Rest: Vgsr). Carbon stars without 
radial velocities are illustrated by red crosses at V(3Sij=333. Middle panel: RA against heliocentric distance. Stars with a heliocentric 
distance > 100 kpc are marked with outer circles - it can be seen that these lie exclusively near the Sgr orbital plane. Bottom panel: 
RA against Dec. The red dotted line is the Galactic Plane and the black dotted lines are ±15°; within this region of the Galactic Plane, 
Catalina Survey data are not available The limits in declination also correspond to the Catalina limits. The approximate regions of 
TriAnd (blue polygon) and Hercules-Aquila Cloud (black polygon) are also shown. The blue dotted curves show the ecliptic, ±10°, 
indicating the region of the MLS survey used bylDrake et al.l (|2013bl. 
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Figure 16. Location of carbon LPVs in relation to the LMIO model (small points) for the last wrap; showing the leading (green) and 
trailing (red) arms, as seen in X, Y and Z coordinates. The Sun is 8.5 kpc from the Galactic Centre. The dotted black lines illustrate the 
cone of avoidance due to the Catalina Surveys reaching down to ±15°. Those stars attributed to key halo substructures, in addition to 
Sgr, are indicated with surrounding symbols. For clarity, we only show the region within 120kpc of the Galactic Centre, hence a few very 
distant stars are not seen on these plots. Two stars at [X,Z]~[-50,-50] are probably C sequence stars that have been wrongly classified 
as C', by our selection cuts. 


Catalina data, at RA from 350-360°. The large extent of A16 
makes it more plausible that onr carbon stars are associated 
with this structure. It is possible that HG118 may also be 
part of A16, but due to its much greater distance (~180 
kpc) and uncertaintly over its nature we do not include it 
as such. 


4-1-3 Gemini Arm (Feature 2) 

Feature 2 comprises two distant carbon stars (HG39 and 
the closer HG45) that sit near a feature first found in BHBs 
by Newberg et al.| ( |2003 1, further studied in|Ruhland et alT 


(20111, and later denoted the “Gemini Arm” by Drake et 
|al. I ( |2013^ . [Belokurov et al.| ( |2014p , argue that the Gemini 
Arm is actually the true location of the apocentre of the 
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Figure 17. Spectra for distant (R >100 kpc) carbon stars with 
SDSS data. The left panels show the full spectrum for 4000 - 
9000A, and the right panels show the corresponding region for 
4000-5000A, where the G-band sits at ~4300A. 


objects is rather difiFuse, and appears to extend from the 
known apocentre of the Sgr arm to a distance of more than 
100 kpc. These stars may be part of the Sgr leading arm 
proper, or they may be a distinct, more distant feature. 

A hint of an extension to the leading arm was seen in 
[Belokurov et ah' (20061, and the other tracers suggest that 
the situation at this location is complex. In the central panel 
of Figure jlSj RR Lyrae from Drake et al. (2013b I can be seen 
at this location. Again - as noted earlier - these are from the 
MLS which has limited sky coverage, and the actual number 
of RR Lyrae at this distance may well be greater. 

Further tracers in this region include two giant stars 
from Dohm-Palmer et al. (20011, which they suggest may 
have come from an older wrap of the Sgr streams. There 
are also BHBs ( Clewley et al.|2005 ), six of which lie in two 
velocity gr oups: four at Vcsit^ 20 kms“^, and two at ~- 
50 kms“^. Clewley et al. (20051 argue this represents stars 


at the turn-around of an elliptical orbit, although as their 
BHBs do not lie on the model of |Law, Johnston, fc Majewski| 
(20051, they do not argue for an association with Sgr. 

Another carbon star, HG57, also lies at a large distance 
at this right ascension, but sits far to the south of the other 
tracers. This star is isolated from any known halo structure, 
and so we assign as “Unclear”. 


4.L5 LMC 


trailing arm of Sgr stream, explaining why that part of the 
LMIO model has not been seen at the distance expected in 
the area north of the Galactic Plane. 

The Gemini Arm is particularly apparent in the RR 


Lyrae of Drake et al. (2013b I. They also comment on several 


of the most distant RR Lyrae, which they think may be 
distinct from the main body of the Gemini Arm. These more 
distant RR Lyrae lie close - on the sky - to HG39, and are 
also at about the same distance. It is significant that these 


Star HG32 is some 7.7° on the sky from the centre of the 
LMC, compared to the tidal radius for the LMC of ~ 17.1° 
(van der Marel et al. 20021. It has a distance of ~50kpc. 


which is also consistent, and this star is most likely a mem¬ 
ber of the LMC. Likewise, at 12.2°, HG35 also a probable 
member of the LMC. HG31 lies 18.1° away from the LMC, 
just beyond the tidal radius, but is 55 kpc distant and may 


also be associated. Munoz et al. (20061 found stars 22° from 


the centre of the LMC, near the Carina dwarf, and HG31 
lies in the direction between the LMC and Carina. It would 
be valuable to obtain the velocity of the carbon stars to 


RR Lyrae come from the MLS component (which reaches 

see they are similar to the giants that Munoz et al. 

(20061 

much deeper than the GSS and SSS surveys) of the Catalina 

proposes to be a part of the LMC. Majewski et al. 

(20091 

Surveys, but which is only 20° wide, straddling the ecliptic. 

found RGB stars out to the limit of their survey - 23 

° from 


These distant RR Lyrae stretch across the full width of the 
MLS strip, and so likely cover a wider region of the sky, 
possibly including the position of HG39. 


Sesar et al. (20121 report two kinematically distinct, but 


spatially close, groups of RR Lyrae stars (denoted Cancer 
A and Cancer B). They both lie at the same location in the 
sky and distance (~86 kpc), but differ in radial velocity: the 
RR Lyrae of Cancer A have a Vgsr of 78 kms“^ (ct = 12 
kms“^), while those in Cancer B have a Vgsr of 16.3 kms“^ 
((7 = 15 kms“^. The Vgsr of HG45 (114T10 kms“'^) might 
argue against a connection with the Gemini Arm, but it may 
be related to Cancer A, given the offset in distance between 
it and the RR Lyrae. 


the centre of the LMC, which they attribute to a “classical” 
halo. 

If the carbon stars are indeed part of the LMC is it 
an interesting result, as it would suggest that either the old 
LMC halo has an intermediate-age component, or that such 
stars have been stripped from the LMC disk by interac¬ 
tions (Munoz et al. 20061 - although the latter scenario is 


challenged by work which indicates the LMC disk is un¬ 
perturbed, at least out to a galactocentric radius of 16 kpc 
( Saha et al.|2010 1. 


4-1.6 The Sgr “bundle’’ 


4-1-4 Sgr apocentre (Feature 3) 

Feature 3 lies near - but somewhat more distant than - the 
apocentre of the leading arm of the Sgr stream. It comprises 
only one carbon star (HG66) at a distance of 125 kpc, but 
is surrounded by a number of other tracers. This group of 


There is a “bundle” of LPV carbon stars, very close to each 
other on the sky, (at RA, Dec, ~ 300°,-30°). These lie on 
the outer region of the main body of the Sgr dSph, although 
the main part of the Sgr galaxy itself is close to the Galactic 
plane, and so not within the Gatalina Surveys footprint. 
These are almost certainly Sgr stars. [McGonnachie (20121 
gives the half-light radius of Sgr as ^^5.7°, while [Majewskij 
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et al. (20031 derive a King fit with a limiting radius of 30°, 


which would include these stars. 


4-1.7 Galactic Anti-centre structures (GAS) 


There are a number of known features in the region of the 
Galactic Anti-Centre, especially north of the Galactic Plane. 


These include the Anti-Centre Stream (ACS) of Grillmair, 
[Carlin, fc Majewski|(|2008| ), and the Monoceros Overdensity 
(Newberg et al. 20021. Of these only Monoceros has been 
detected south of the Plane (Meisner et al.||2012 Slater et 

[S)][M4| ). 

Two of our stars (HG40, HG41) stand out as having 
velocities inconsistent with the LMIO model. Both are near 
the Galactic Plane, are relatively close to the Sun (11, 14 
kpc); and have locations, distances and velocities consistent 
with of the Anti-Centre Stream. Another, HG36, although 
it lies close on the sky to the trailing arm of Sgr, is too close 
at 12 kpc to be part of the stream. All three are consistent 


with of the Anti Centre Stream as found in Li et al. (20121 


One carbon star in the south, HG29, of the Galactic 
Plane may belong to the Monoceros Overdensity. It just lies 
in the region discussed by ( Slater et al.|2014 |, at a distance 
of ~12 kpc. 


4.1.8 TriAnd 


At a location of approximately RA ~30° and Dec ~35°, we 
find a small group of carbon stars (HGl, HG7, HG8, HG12, 
HG14 and HG15) with a distance of about ~20kpc. This 
region corresponds to the Triangulum-Andromeda (TriAnd) 
structure, a relatively dispersed stellar system found in M 
giants ([Rocha-Pinto et al.|2004|) , but which has also been de 




tected in MS and MSTO stars ( [Majewski et al. 2004 Martin, 
Ibata, fc Irwin|2007 ( (blue polygon on lower panel of Figure 


The TriAnd feature, as shown in Figure 2 of [Rocha-[ 


Pinto et al. (20041, has a suggestion of two separate regions. 


and five of the six candidate TriAnd carbon stars all lie in 
one of these regions, only HGl in the second. This may indi¬ 
cate that the two regions have different stellar populations, 
although the small number of carbon stars involved prevents 
any firm conclusion. The location of our carbon stars to one 
region of TriAnd is not due to a lack of Catalina Surveys 
data at low right ascension. 


Sheffield et al. (20141 note that Martin, Ibata, & Irwin 


(20071 report two features, denoted TriAndl and TriAnd2, 


which they find to have distances of 15-21 and 24-32 kpc 
respectively, with ages of 6-10 and 10-12 Gyr . The first 
sample carbon stars are more consistent with TriAndl, from 
distance, and this is also consistent with the expectations of 
an intermediate-age population. Only HGl of these carbon 
stars has a published velocity (Vgsr), 69 kms“^, consistent 
with membership of TriAndl as determined bylSheflield et| 

31(|M4l. 


4-1.9 Hercules-Aquila Gloud (HAG) 

The Hercules-Aquila Cloud (HAG), discovered by|Belokurov[ 
[eTS;] ( [2007| ), is a large and diffuse feature straddling the 
Galactic plane (black polygon on lower panel of Figure [l^. 


It stretches ~80°, centred on 1=40°, appears to extend either 
side of the plane by some 50°, and lies at a distance of 10-20 
kpc. In the north, Belokurov et al. (20071 have evidence of 
a Vgsr of roughly centred on 180 kms“^. 

RR Lyrae have been found that correspond to the HAG 
by both Watkins et al. (20091 and Sesar et al. (2010a I in 
SDSS stripe 82. More recently, [Simion et al.[ ( [2014[ ), also us¬ 
ing RR Lyrae as a tracer, find a strong asymmetry in the 
HAG. It is very prominent in the Galactic Southern Hemi¬ 
sphere, but barely seen in the North. But they do suggest 
that this is due to interstellar dust reaching to greater Galac¬ 
tic latitudes in the Norther Hemisphere. 

HG91,HG102, HG104, HGl 10 and HGl 12 are good can¬ 
didates for this feature, especially we consider a region of the 
sky somewhat larger than drawn on Figure [T^ Two others, 
HGHG80 and HG85 might be associated, and we tentatively 
allocate them to the HAG. 


4-1.10 Starkenburg Group 6 

HG52 (RA,Dec ~ 170°,-18°) lies off the LMIO model - as 
seen on the sky - far to its south. This star has a distance 
of ~40 kpc and a velocity of 169 kms“^. It is 
of the pairs of K-giants (their group 6) from 

on the sky, and has a similar distance and 
short-dashed box in all panels of Figure [T^ . 
The Starkenburg pairs are good candidates to be tracing 
halo substructure due to the improbability of random halo 
stars sharing location, distance and velocity. The presence 
of HG52 supports the reality of the group. 

HG48 may also be associated with the Starkenburg 
group 6, as it too has an unusual radial velocity that places 
it close to both this group and to HG52. 


et al. (2009 


velocity (see 


close to one 
Starkenburg 


4.2 Notes on specific stars 

During our analysis, a number of individual stars stood out 
as worthy of comment: 


4-2.1 HG77 - SDSS J144945.37-h012656.2 


One star [HG77] is unusual in that it has a strong GH feature 
in its spectrum (Figure [l^ , is quite blue, with a (J — Ka)o 
colour of 0.9 mag, and a short (88 day) period. Although 
with a GH feature, this star is probably a G-R star. Gom- 
parison with the criteria listed by |GoswarnI| ( [2005[ ) - the 
sharper Na D feature and strong GN lines - are indicative 
of a G-R star, as is the weakness of Hq . 

Our belief that this is a late-R type AGB, rather than 
a GH-type star, is strengthened by its location (assuming 
it is an AGB star). It position on the sky, and its derived 
distance (52 kpc) place it in the apocentre of the Sgr leading 
arm, accompanied by many other carbon stars. 


4-2.2 Star S25 

The star S25, found in Sext ans ([Mauron et al.[[2004p , has 
been previously reported by Sakamoto et al. (20121, using 
their own photometry, which has fewer data points than in 
the Gatalina Surveys. We analysed this star in more detail, 
as a carbon star with such a long period (and hence not old) 
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Figure 18. Spectrum for star HG77, from SDSS data. The inset 
shows the blue region of the spectrum, and the CH feature at 
~4300A. 



Figure 19. Plot of all LINEAR and CSS data for the star S25, 
with a model curve that is the best fit for all the Catalina data 
(red). The red vertical dashed lines bound the obscuration event. 
We also show the I-band data from [Sakamoto et al!] )|2012|l . 


is not expected in this dwarf galaxy, which is metal-poor and 


primarily composed of an ancient stellar population (Lee 


et al. 20091. The light curve and initial fit of S25 are also 


shown in the Appendix, to be consistent with the full sample. 
However, it was sufficiently of interest to warrant further 
analysis. 

The extensive Catalina data have been supplemented 
LINEAR observations, to reveal some of the behaviour prior 
to the Catalina data being taken. The LINEAR and CSS 
data overlap between MJD = 53800 and 54600, and within 
the range MJD= 54400 and 54600 they both have small er¬ 
rors in both samples. We use this region to determine the 
offset between the two systems as 0.1 magnitudes. This cor¬ 
rection is then applied to the CSS data for the subsequent 
analysis. 

The light curve (Figure |19[ ) shows it suffers from an 
obscuration event in the the region from MJD 54,600-55,900 
days (between the red vertical lines), with a drop of over four 



Figure 20. Light curve for HG2-42 with the best Period04 model 
fit, with a period of 1744 days and an amplitude of 1.8 mag. 


magnitudes. A fit to all the Catalina Surveys data (the red 
curve in Figure [l^ seems to explain the obscuration event as 
actually being part of a longer secondary period. However, 
we can see that this fit does not match the LINEAR data 
well. Thus we believe that this is a genuine obscuration event 
from significant mass-loss. 

Even we exclude the ’dip’, with its larger uncertain¬ 
ties, the total range in amplitude of the combined LINEAR 
and Catalina Surveys data is still over 3.8 magnitudes. This 
compares to a range of ~3 mags for SI, (M29 in the Mau- 
ron numbering system). Sloan et al. ( 2009| ) found that M29 
in the Sculptor dSph is undergoing substantial mass-loss, 
and if S25 is similar it will provide an additional example 
of the production of carbonaceous dust in a metal-poor en¬ 
vironment. Hence, S25 is worthy of observational follow-up. 
Additional photometry will add to the light-curve data to 
better understand the nature of the long-term changes in 
luminosity, i.e. whether they are periodic or not. 


4.2.3 HG2-42: [TI98]1450-1300 

HG2-42, originally from [Totten fc Irwin (19981, has a very 
long period of 1638 days and an amplitude of 1.8 magnitudes 
(See Figure |20| ). We exclude from our first sample, as this 
period lies far from the region in which we could expect the 
period-luminosity relation to hold. Hence it is not in Table 
but we place it in Table using its luminosity to estimate 
a distance. 

It has been suggested that stars having periods longer 
than ~1000 days are massive enough for hot bottom burning 
(HBB) to occur ( Whitelock et al)]|2013 l. This will result in 
carbon being consumed and the carbon star turning into an 
O-rich star ([Feast|[7009|). The carbon stars with the longest 


period in the LMC are just below 1000 days (Whitelock et 
al.|2003 1. So this carbon star is quite unexpected. 


Nanni et al. (20131 have models in which a very low 


metallicity star may become C-rich if it is massive (M > 
5Mq), due to a very efficient HBB in which the ON cycle 
is activated, leading to the depletion of oxygen ( SiesspoTo |. 
However, the Vgsr of HG2-42 places it in the Sgr leading 
arm, and the distance derived using NIR photometry, is also 
consistent with a location in the Sgr leading arm. Although 
circumstantial, this does suggest that it may not be a very 
massive star, but a normal AGB star with a very long period. 

This idea may find some support from jSoszynski et al.j 
( [2009[ ) who find a group of 10 Miras with P>1000 and am¬ 
plitudes in the I-band of 4mag. As they exhibit ’regular’ 
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variability, they assume they are 0-rich, not C-rich. How¬ 
ever, we found that is not the case for shorter period stars 
(in that many of our carbon stars show very regular light 
curves), so maybe there are carbon stars among this group 
of 10. By extension there may be a population of very long 
period carbon Miras that have hitherto not be found due to 
their scarcity. 


5 CARBON STARS WITHOUT PERIODS 


For many stars in our input lists (see we do not have pe¬ 
riods. This can occur for a number of reasons. In many cases 
the star lies outside the Catalina Surveys footprint, in others 
the Catalina light curves have too few (or too noisy) data- 
points to obtain a period fit, or show no obvious variability. 
And finally, in some cases, the light curves exhibit behaviour 
that prevents a period fit being made. This is often the case 
for very bright stars in which the Catalina data are satu¬ 
rated. For these stars we use the more common method of 
distance estimation from NIR photometry alone. 

We were encouraged to do so when we found that the 
distances obtained by Mauron et al. using Ag-band magni¬ 
tudes and {J — Ka) colours were relatively consistent with 
those derived by using periods. This consistency may, of 
course, only be the case for those AGB carbon stars which 
have well-behaved light curves, or at least are shown to be 
variable. Therefore although we include the stars without 
periods in this section of the paper, we caution against 
drawing very strong conclusions from them. We include a 
list of them here, to encourage further observations of those 
with weak light-curves, or none at all (i.e. those not in the 
Catalina footprint) 


We remove stars classihed as CH-type by 

Totten & Ir- 

win ( 19981), and those in the list of CH-stars of 

Bartkevicius 


(19961, as these will not fit the relation that is true of AGB 


stars. Of course in some of our sources, an object is reported 
as a carbon star, without further details as to the sub-type, 
so one must be further cautious of the results from the sec¬ 
ond sample. To address this problem, we also exclude those 
stars for which (J — Ks)q < 1.2, as such blue colours are 
typical of CH-stars (e.g. see Fig.[4|), unless they are noted as 


C(N) by Gigoyan fc Mickaelian| (20121. Finally, we remove 
stars for which Ks < 6, as these will all be in the disc. This 
new sample, which we term the “second sample” is shown 
in Tabled 

We now estimate the distances of the second sample 
stars following the approach of Mauron et al. They exploit 
relations between absolute Ks -band magnitude and J and 
K photometry, derived from a sample of LMC carbon stars. 
The equation(s) used by Mauron et al. to estimate distances 
are not given, so we re-derive them from the data given in 
their papers. We use the reddening corrected NIR photome¬ 
try and derived absolute Ag-band magnitudes from |Mauron,| 
Gigoyan, & Kendall (2007b \ and Mauron (20081 to obtain 


two fits: 


Mks = 1-550 X (J - As)o - 4.927 for 1.2 < (J - Ks)o < 1.8 

( 11 ) 

Mks = 0.414 X (J - Ks)o - 8.429 for (J - As)o > 1.8 

( 12 ) 



Distance (kpc) - from NIR 


Figure 21. Plot of our distances, derived from variable periods, 
against their distance determination from NIR photometry. There 
is a reasonable correlation between the two distance measures for 
this sub-sample. The line of equality (dotted) is shown to guide 
the eye. 


We can test these in more detail by comparing the dis¬ 
tances obtained from periods, against those from NIR pho¬ 
tometry alone, for our first sample (see Fig. 211. Given the 
distance errors the correlation is remarkably good. A few 
outliers are outside of the estimated distance uncertainties: 
HG118, HG77 and HG106. The hrst two of these have given 
us problems before: HG118 has a weak light curve and was 
found to be an outlier in Fig. |14[ HG77 is one of the bluer 
stars with a very short period. HG106 also has a weak light 
curve. 

Equations and |12| are then used to derive estimates 
for Mk, and then distances for the stars in our second sam¬ 
ple. This approach effectively assumes that the carbon stars 
are all on the C sequence. Hence any that are, in fact, C' 
or D will have incorrect distances. But without period and 
amplitude data, it is difficult to know which sequence our 
sample stars belong to. The exception to this occurs if the 
(J — Ag) colour is greater ~2, in which case it is reasonable 
to assume that it lies on the C sequence (see Figure]^. Some 
stars in this sample have published velocities, which are also 
given in Table and can increase our conhdence of whether 
they lie in known halo substructures. 

The good correlation in the two methods of distance 
estimation seen in Fig. |21|may, of course, be a result of the 
nature of sub-sample of carbon stars that are found in the 
hrst sample. That is, the hrst sample have variability prop¬ 
erties that make the use of NIR photometry appropriate. 
Many stars in the second sample, however, show either no, 
or erratic, variability. For this reason, we repeat our note of 
caution from above. Although we believe that the distances 
are indicative for the sample as a whole, individual stars 
may have estimated distances that are incorrect. Due the 
additional uncertainties of this sample, we do not attribute 
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the second sample stars to substructure as done for the hrst 
(LPV) sample. 

Looking at the position of this second sample in ve¬ 
locity, distance and location (Figure \22\ a few results are 
immediately apparent. Firstly, many more stars (shown as 
asterisk symbols in the figure) than our “first” sample have 
positions that place them close to the MW disk. This is be¬ 
cause many of these stars are too luminous to obtain good 
light-curves in the Catalina Surveys; they become saturated, 
making curve htting impossible. 

We also see many stars that are almost certainly in the 
LMC (either side of RA, Dec ~ 90°, -70°) and SMC (sur¬ 
rounding RA, Dec ~ 15°, -70°). Most of these have no light- 
curves as they are in the region excluded from the Catalina 
Surveys due to crowding in the imaging. Others are beyond 
the southern declination limit of the Surveys. The second 
sample carbon stars nearby to the SMC reach to a projected 
sky distance of 4.3°, equivalent to ~4.7 kpc, a value consis¬ 
tent with the recent results of [Noel fc Gallart] ( |2007| ), who 
found evidence of an intermediate-age population up to 6.5 
kpc from the centre of the SMC. So these stars are consis¬ 
tent. 

One interesting feature are the several “second sample” 
carbon stars at distances >100 kpc, and that lie close to 
the features identified in the first sample. One lies close to 
the A16/Pisces Overdensity, at a distance of over 150 kpc, in 
“Feature 1”. Three second sample stars lie beyond the apoc- 
entre of the Sgr leading arm, from 100-160 kpc in “Feature 
3”; two on the LMIO model (on the sky), and the third is 
not far from it. 


merger activity of massive dwarf galaxies onto the MW, ex¬ 
cept for the Sgr dwarf, and any progenitors of the features 
discussed above. Of course, our LPV sample in particular, 
avoids the Galactic Plane, so any accretion events occurring 
in the Plane would not be detected. 

This hnding is consistent with the view that the MW 


has had a quiescent history (Hammer et al. 20071. Specih- 


cally, accretion into the field halo by satellites with extended 
SFH has been unimportant for the last ~8 Gyr (Wyse 20081. 


However, older, more metal-poor satellites may well have 
been accreted, or are still being so, in the more distant halo. 

However, we find - as many others have using a variety 
of tracers - that the LMIO model is not a good ht to the 
data along its full spatial extent. There does indeed appear 
to be a lack of carbon stars in the trailing arm north of 
the Galactic Plane, where the model suggest a relatively 
high stellar density (and as noted above). Only HG2-26 is 
plausible as being consistent with the LMIO, but - as noted 
earlier - it could equally be part of the Gemini Arm. The 
lack of the trailing arm in the region predicted by the LMIO 
model has been reported before. For example, [Drake et al.| 
1 2013a I have also found that their RR Lyrae show little of 
the LMIO trailing arm at RA = 110° at a distance of ~40 
kpc. 

But our data are also ambiguous with regards to the 
leading arm south of the Plane. Although we have carbon 
stars that we attribute to Sgr in this region, they are few. 
Moreover, many of these may in fact be part of the Mono- 
ceros Overdensity. Such a result, as with the trailing arm in 
the north, is not totally new. Other tracers have either not 


Finally, HG2-26 is particularly interesting. Its location 

been detected in the leading arm (e.g. 

Sesar, Juric, & Ivezic 

on the sky (RA, Dec ~ 106°, 24°), and its estimated dis- 

2011b||Ruhland et al.||2011||Sirko et al.||2004||Carrell, Wil- 

tance of 65 kpc would place it in the region of the trailing 

helm, & Chen|2012l, or have been only detected to a minor 

arm in which no other tracers have been found. It is in the 

extent (Shi et al.||2012 1, in this region of the sky. 


Catalina data, but with too few data points to obtain a good 
light curve. If we use the same distance uncertainty used for 
the hrst sample, based on the (J — K^) colour, this star 
could be at a distance from ~47-90 kpc. The lower value 
would be consistent with the LMIO model, while the upper 
value would locate it in the Gemini Arm (which of course 
may be the true Sgr trailing arm). The Catalina data for 
this star range over 3 magnitudes in the V-band, suggest¬ 
ing it also suffers from episodes of obscuration, and that the 
uncertainty in the distance may be greater still. 


6 DISCUSSION 

In the following discussion we focus on the results from the 
hrst sample (Table [^, rather than those for which distances 
have only been obtained from magnitude (Table [^, due to 
the far greater uncertainties associated with second sample. 


6.1 Sgr carbon stars 

A major hnding of this study is that the majority of the 
hrst (LPV) sample can be understood as being part of the 
Sgr streams, being consistent in location, distance and ve¬ 
locity (where available). This highlights the importance of 
the Sgr accretion in the history of the MW halo, especially 
as seen through intermediate-age tracers, such as carbon 
stars. This result suggests that there has been little recent 


6.2 Distant carbon stars 

Another result is the presence of a handful of carbon stars 
at very large distances, greater than 100 kpc. Feature 1 is 
almost certainly part of the Pisces Overdensity, and Feature 
2 is very likely tied up with the “Gemini Arm”. Feature 3 is 
more difficult to explain. As these stars lie on the Sgr orbital 
plane, one plausible explanation is that they arise from older 
wraps of the Sgr streams, lost on previous orbits, especially 
if the Gemini Arm is indeed the actual apocentre of the Sgr 
trailing arm. 

The idea that the most distant carbon stars may be 


older wraps of Sgr is consistent with Shi et al. (20121 who 


argues that it lost an older, more metal-poor, component in 
the older wraps, and then the younger more metal-rich com¬ 
ponent in the youngest wraps (the arms shown in Figures [T^ 
and \22\ . For example, HG39 has a period of 210 days and 
HG66 of 149 days. The shorter periods are indicative of a 
greater age than the majority of the sample. 

An alternative scenario ( [Garlin et al.|[?012p suggests 
that both the Pisces Overdensity and the Gemini Arm can 
be best explained as the tidal features of the progenitor of 
the Virgo Overdensity; a feature that is located at a distance 
of ~ 15 kpc, over a large area of sky of the SDSS footprint 
above the Galactic Plane. Guriously, we have no obvious 
candidate carbon stars for this structure, whereas we have 
examples in both the Pisces Overdensity and Gemini Arm. 
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Figure 22. A similar plot as Figure [T^ but including the second sample (with no periods). The symbols indicating possible feature 
membership of the first sample are the same as those shown in Figure [T^ The RR Lyrae data of [Drake et al.| l |2013b| l are excluded to 
for clarity. 
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An obvious issue concerns the reliability of the distance 
determination for the handful of very distant carbon stars. 
One concern might be that they lie in the plane of the Sgr 
orbit because they may actually be in the main arms of the 
streams (i.e. up to 60 kpc) and are just projected to larger 
distances by failings in our approach. However, these stars 
have well-behaved light curves and bluish colours. As dis¬ 
cussed above, the spectra of four of the most distant stars are 
N-type. We doubt that the NIR data have been affected by 
obscuration events that may confound the P-L relationships; 
there is no evidence that these suffer erratic and dramatic 
changes in their light-curves (as seen in S25), in which the 
2MASS NIR data would not be representative of the period 
covered by the Catalina Surveys light curves. The long dura¬ 
tion of the Catalina Surveys, some seven years, is invaluable 
in assessing this. 

By contrast, as we have seen above, other carbon stars 
were difficult to assign cleanly to the C or C' sequences, and 
some proportion of them may be closer, if they are indeed 
C' stars. However, as we are interested in using these stars 
as potential tracers of distant structure, this issue can be 
addressed through follow-up observations. Deep imaging in 
the regions near these objects should reveal any halo features 
at the indicated distances. 


6.3 Other Substructure 

In addition to the many carbon stars in Sgr, and the distant 
stars which are in halo structures which may or may also not 
be associated with Sgr, we found LPV carbon stars in other 
known substructures. In particular, we find a clear group in 
Triangulum-Andromeda feature, and several that are likely 
in the GAS. Less clear are the stars we attribute to the HAG, 
primarily as the true extent of this feature is still uncertain. 

The most interesting feature, for our purposes, are the 
two carbon stars near to the Starkenburg group 6. Even 
though we assign these carbon stars as “Unclear”, the group¬ 
ing in velocity, distance and location strongly suggest that 
this is a genuine halo feature. Follow-up imaging in this re¬ 
gion should provide a definitive answer. This is the best 
candidate for “new” substructure. 

Finally, there are several isolated carbon stars which 
are not close to any known substructure. These may be con¬ 
taminants (e.g. non AGB stars; possibly HG107, for exam¬ 
ple) that have fallen into our sample due to our generous 
selection approach, or they may indicate unknown halo sub¬ 
structure. What is remarkable is that there are so few. This 
is even more so when one considers that some of the stars 
that are not easily attributed to a substructure are also those 
with the most uncertain light curves (e.g. HG64 and HG72). 


6.4 Period distribution of tbe carbon stars 


One of our conclusions, that the majority of the carbon stars 
are associated with Sgr, is quite different from that of|Bat-| 


tinelli & Demers (20121. They argue that the halo carbon 


stars were too young to be from the Sgr streams, based on 
the distribution of periods for their sample. They argued 
that their sample has too many ’shorter’ period LPVs than 
would be expected if the Sgr dSph were a relatively old to 
intermediate-age system. 
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Figure 23. Histogram of periods for all our sample (hatched), 
those attributed to Sagittarius (labelled ’Sgr’ and ’Sgr:’ in the 
tables; cross-hatched), and those with a declination less than - 
20° (solid fill), comparable to the sample of |Battinelli fc Demers] 
l|2012[l. The vertical dashed lines delineate the ages adopted by 
[Battinelli fc Demers] ([2012^. 


But we suggest that their conclusion is based on a spa¬ 
tially biased sample. In Figure we plot the distribution 
of periods for our first sample (single hatched), and those 
we attribute to Sgr (cross-hatched). The peak for the Sgr 
stars is at 200 - 250 days, not the 250 - 350 days found by 


Battinelli & Demers (20121. However, they had an implicit 


selection, due to observational constraints, of a declination 
less than -20°. If we apply this same selection to the first 
sample (solid fill), we too find a shift in the peak to the longer 
periods, 300-350 days, consistent with |Battinelli fc Demers] 
(20121. This shift in period distribution is most likely due 
to the large proportion in this sample of carbon stars that 
lie close the main body of Sgr. There is known to be an age- 
metallicity gradient along the streams ( |Ghou et al.|2007]|Shi| 


et al. 2012), so the stars near Sgr have been recently stripped 


from the main body and so represent a younger population 
than the bulk of the Sgr stream stars. 

However, the age of the carbon stars is unclear. The 
short periods of many of the carbon stars should indicate 
an older age. A plot of Mk^ against (J — Ks) colour (Fig¬ 
ure [2^ shows where the LPV samples sit in relation to the 
isochrones of Marigo et al. ( 2008[ p The majority of the sam¬ 
ple lies between the isochrones (with a metallicity of [Fe/H] 
= - 0.7, an average value for the Sgr body and arms, [Oiou| 


et al. 20071 for 1.5 and 2.0 Gyr, suggesting that the sample is 


relatively young. Alternatively, we could be seeing the conse¬ 
quences of the wide range in metallicity (~ 1 dex) observed 
in Sgr and its arms. The total sample looks very similar in 
distribution to the carbon stars in the Fornax dSph, which 


® http://stev.oapd.inaf.it/cgi-bin/cmd 
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Figure 24. Plot of absolute it’s-band magnitude against (J-K 5 ) 
colour. Red squares indicate stars that lie in our ’disklike’ re¬ 
gion (see text), and green squares indicate stars that are in Sgr. 
Also shown are isochrones from the tool CMD 2.5 ( [Marigo et al.| 
|2008[ |, for a metallicity (Z) of 0.004 ([Fe/H] = -0.7) and ages of 
1.5 (green) and 2.0 (blue) Gyr, using the bolometric corrections 
of |Loidl, Langon, &c j 0 rgensen| ( |2001| and dust compositions of 
0.85 AMC (amorphous carbon) + 0.15 SiC for C-rich stars and 
0.6 AlOx-l-0.4 silicate for the O-rich stars ( |Groenewegen|[2b06[ | . 
The horizontal dashed lines shows the tip of the red giant branch 
(TRGB) value of for [Fe/H] = 0.0 (upper) and [Fe/H] = —2.5 
(lower) from [Ferraro et ^ l|2000||. The outlier HG107 is labelled. 


had a fairly continuous star formation, from 14 Gyr ago un¬ 
til as recently as 0.25 Gyr (Grebel fc Stetson|jl999 de Boer 


jet al.|2012[ ). It is possible that the process of accretion trig¬ 
gered star-formation in the Sgr dSph over the last few Gyr 
created both a spread in metallicities and ages, or that it 
contained a range of populations at the time of accretion. 

It is not obvious why this inconsistency between the 
ages implied by the period of variability and the colour- 
magnitude diagram occurs, although there are still unsolved 
issues in AGB evolution, and the models may require further 
refinement. Recently Girardi et al. (20131 found that, from 
their models, the numbers of TP-AGB stars are “boosted” 
by a factor of ~2 for ages from 1.57 to 1.66 Gyr. This may 
help explain the the presence of so many of our sample (and 
the MW satellite carbon stars) seen in Figure [24{ and which 
sit between the 1.5 and 2.0 Gyr isochrones. 


7 SUMMARY 

We have revisited the use of carbon stars as tracers of halo 
substructure, with a particular focus on the Sagittarius dSph 
tidal streams. 

Garbon stars on the AGB have always promised the 
opportunity to trace intermediate-age stellar populations to 


large distances in the MW halcj^ Any such features will 
likely be the remnants of disrupted satellite galaxies. The 
potential of these stars has been difficult to realise as they 
can have a wide range of luminosities, making distance de¬ 
termination imprecise. By using variable carbon stars, we 
exploit period-luminosity relations to obtain new estimates 
of their distances. Our own sample suffers from the inherent 
uncertainty in our NIR photometry, due to it being single 
epoch only, but in many cases these distance uncertainties 
are comparable to those found in the use of other tracer 
populations. 

This study been made possible by the availability of 
the recently released Catalina Surveys Data Release 2 light 
curves. These cover a large part of the sky and have data 
over a duration of seven years, allowing analysis of stars 
whose periods can be more than a year. The relatively large 
number of carbon stars that we have assembled from the lit¬ 
erature means that those with inadequate light-curves can 
be excluded from our “first sample”, such that we have con- 
hdence in the results from this sample of LPV carbon stars. 

Almost all of the stars in our “first sample” can be 
allocated to known substructures. The majority can be rea¬ 
sonably assigned to the tidal streams of Sagittarius. Other 
features that possess AGB carbon stars are the Hercules- 
Aquila Cloud, the Galactic Anti-centre streams, the LMC 
and Triangulum-Andromeda. 

We also find that the LMIO model is consistent with 
many of our Sgr carbon stars. However, we also find two 
possible regions of failure in the model. There is a distinct 
lack of carbon stars in the trailing arm north of the Galac¬ 
tic Plane (except for the highly-schrouded HG2-26 and its 
associated large uncertainty in estimated distance), and ev¬ 
idence of the leading arm south of the Plane is also weak. 
Due to their location near to the Galactic plane (with the 
consequent crowding and high extinction), it is difficult to 
draw very strong conclusions, but these apparent carbon 
star dehciencies remain a challenge for the model. The re¬ 
gions involved are also confused by the presence of other 
known substructures: the Anti-Centre Streams, Monoceros, 
and the Hercules-Aquila Cloud. Future observations to ob¬ 
tain velocities for many of these carbon stars may be valu¬ 
able in identifying their parent. As crucial will be data that 
can probe closer to the Galactic Plane than hitherto, as it is 
in this region that the fate of both the leading and trailing 
arms should become clearer. 

Despite two decades of study, the Sgr accretion still 
poses many questions. The bifurcation is still unexplained, 
and the Gemini Arm has been suggested as related to ei¬ 
ther the Sgr dSph or to the Virgo Overdensity. It is no sur¬ 
prise that Koposov, Belokurov, & Wyn Evans (20131 remark 
that the “ Sagittarius (Sgr) dwarf galaxy remains a riddle, 
wrapped in a mystery, inside an enigma”. 

Carbon stars provide new probes in the far halo which 
may point to solutions to these problems. They are too few 
in number to say much on their own. But if the distant halo 


® However, the need for caution is highlighted by the recent find¬ 
ings of jFeast, Menzies, &: Whitelock] | |2013| l who found a carbon 
Mira with a period of 551 days in the globular cluster Lynga 7. 
Despite its long period, it must be old, and they suggest it is the 
result of a stellar merger. 
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carbon stars are representatives of a wider stellar population 
then deep imaging of the regions around these stars may re¬ 
veal new stellar substructure, and its nature if present. In 
particular, imaging of the regions around the very distant 
carbon stars, and those few that lie away from known sub¬ 
structure would be invaluable. 

Notably, we find that the Trianglum-Andromeda fea¬ 
ture (TriAnd) stands out prominently in our data, and we 
strengthen the likelihood that the Starkenburg group 6 in¬ 
dicates a bone fide halo feature. Very few of our hrst sample 
stars are located in regions where no known substructure is 
thought to exist. One aim of this paper is to identify pos¬ 
sible targets for follow-up, hence we have taken a generous 
approach to the AGB star classification, as they add little 
overhead to follow-up, yet may indicate substructure which 
will show up in deep imaging. The small number of stars 
with no attribution facilitates this. 

We have also highlighted several stars which possess un¬ 
usual properties. Further observations of these objects (e.g. 
photometric monitoring or high-resolution spectroscopy) 
might throw new light on them, and thus on the carbon 
star population more generally. 

In conclusion, we suspect that almost all of major 
intermediate-age substructure in the Milky Way have been 
found already. There may, of course, be features hiding be¬ 
hind the region of the Galactic Plane or towards the Galactic 
Poles, which are not properly covered by our study. The car¬ 
bon star population has also shown that it can still surprise, 
with the complex nature of these objects still able to provide 
examples that ht no easy classihcation. But the the luminos¬ 
ity of the AGB population still makes them one of the best 
probes of the MW halo. 
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Table 1. Catalogue of variable carbon stars. The coordinates are given in decimal degrees (J2000), the (reddening corrected) 
NIR magnitudes and colour are from 2MASS data. The periods and amplitudes (in Catalina V band) of variability are 
derived from fits using Period04. 
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56l? 


222 

0.60 

c 

20 

LMC 

HG32 

103.6610 

-71.5165 

10.95 

2.06 

-7.55 

51+17 


282 

0.89 

c 

6 

LMC 

HG33 

105.4540 

39.9977 

7.58 

2.75 

-7.55 

11+5 


354 

1.22 

c 

14 

Unc 

HG34 

107.7400 

47.9719 

7.33 

1.55 

-7.51 

9ll 


233 

0.47 

c 

18 

Sgr 

HG35 

108.0770 

-63.6360 

11.04 

2.48 

-7.29 

47111 


285 

1.64 

c 

20 

LMC 

HG36 

109.4010 

50.1782 

7.71 

1.65 

-7.73 

1211 

38 

267 

1.29 

c 

1 

GAS: 

HG37 

113.1360 

26.7875 

8.14 

1.71 

-6.50 

9li 


150 

1.26 

c 

14 

Disk 

HG38 

113.6000 

27.3198 

7.53 

1.72 

-7.97 

13l| 


306 

0.79 

c 

3 

Sgr: 

HG39 

119.7330 

33.9772 

13.19 

1.37 

-7.43 

134112 


210 

1.20 

c 

2 

Gemini: 

HG40 

127.3130 

18.3853 

7.06 

1.67 

-8.12 

nil 

18 

323 

1.35 

c 

17 

GAS: 

HG41 

127.3710 

10.7734 

8.13 

2.11 

-7.62 

14l| 

-30 

294 

1.00 

c 

17 

GAS: 

HG42 

130.9810 

-12.4019 

8.02 

1.48 

-7.30 

1211 


204 

0.40 

c 

6 

Disk 
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Table 1. Continued 


ID 

RA 

Dec 

Kso 

(J-Ks)o 

MKso 

R© 


Period 

Amp 

Seq 

Ref 

Feature 


degs 

degs 




kpc 

km/s 

days 

mags 




HG43 

133.5780 

-12.0151 

8.01 

2.71 

-7.79 

+ 1 

-58 

389 

1.13 

C 

17 

Unc 

HG44 

134.3580 

17.3477 

10.70 

4.40 

-6.66 

30t?| 

-105 

394 

1.20 

C 

4 

Sgr 

HG45 

136.4430 

20.4105 

12.28 

2.81 

-6.86 

68ti 

114 

256 

1.58 

C 

2 

Gemini: 

HG46 

138.3830 

19.5729 

7.35 

1.64 

-7.68 

lOti 

-135 

260 

1.25 

C 

16 

Sgr 

HG47 

138.7720 

19.2938 

8.49 

1.68 

-7.65 

17t| 

-104 

258 

1.63 

C 

16 

Sgr 

HG48 

152.6540 

-6.8538 

8.52 

1.41 

-8.37 

24t| 

139 

335 

1.35 

C 

17 

Unc 

HG49 

160.0140 

26.0206 

9.81 

1.32 

-8.11 

38j:f 

-148 

137 

0.22 

C’ 

1 

Sgr 

HG50 

164.8490 

39.7349 

9.45 

1.60 

-7.16 

21±| 

-155 

198 

0.84 

C 

16 

Sgr 

HG51 

167.4990 

-21.3671 

6.68 

1.57 

-8.16 

9tl 

-60 

319 

1.54 

C 

16 

Unc 

HG52 

169.3290 

-17.4876 

11.21 

1.42 

-6.91 

+ 1 

169 

166 

0.86 

C 

16 

Unc 

HG53 

171.7730 

10.2612 

11.66 

1.46 

-6.30 

39t| 


125 

0.64 

C 

19 

Sgr 

HG54 

182.3540 

15.2717 

9.82 

1.33 

-7.87 

35^5 

-81 

256 

0.73 

C 

16 

Sgr 

HG55 

183.5710 

-9.0139 

11.86 

1.33 

-6.86 

56j:f 

-21 

157 

1.04 

C 

17 

Sgr: 

HG56 

186.9170 

-0.4642 

10.54 

2.20 

-7.47 

40115 

-130 

281 

1.76 

C 

1 

Sgr 

HG57 

189.6220 

-40.9361 

12.69 

2.59 

-7.04 

90l|® 


264 

1.83 

C 

15 

Unc 

HG58 

192.2700 

13.3432 

11.16 

1.42 

-7.71 

60l| 

-62 

243 

0.89 

C 

16 

Sgr 

HG59 

192.9580 

1.5005 

10.19 

1.86 

-7.67 

38tV 

-25 

277 

1.39 

C 

1 

Sgr 

HG60 

194.6400 

16.6701 

7.81 

3.47 

-7.30 

111! 

58 

395 

2.00 

C 

1 

Sgr: 

HG61 

195.4830 

8.6087 

12.06 

2.63 

-7.04 

66li^ 


263 

1.21 

C 

2 

Sgr 

HG62 

197.4730 

23.4279 

9.48 

2.26 

-8.65 


62 

505 

1.69 

C 

17 

Sgr 

HG63 

199.7550 

4.3392 

11.84 

1.21 

-6.71 

5ll| 

-2 

141 

0.20 

c 

17 

Sgr 

HG64 

203.3310 

41.9142 

8.17 

1.21 

-7.65 

1511 


221 

0.74 

c 

14 

Unc 

HG65 

203.9880 

6.3986 

10.57 

1.68 

-7.52 

42 If 


242 

1.46 

c 

19 

Sgr 

HG66 

204.6730 

15.8689 

13.45 

0.93 

-7.03 

125llf 


149 

0.61 

c 

2 

Sgr: 

HG67 

206.8460 

-34.7898 

10.08 

1.95 

-7.82 

38lf 

3 

308 

0.85 

c 

17 

Sgr: 

HG68 

208.2560 

0.7874 

11.07 

1.55 

-7.54 

5311? 

-1 

236 

0.64 

c 

1 

Sgr 

HG69 

209.0100 

-1.6073 

11.29 

1.57 

-7.34 

54l| 

-2 

215 

0.68 

c 

16 

Sgr 

HG70 

209.8360 

-30.3940 

11.81 

2.77 

-6.21 

40ll? 

25 

186 

0.92 

c 

16 

Sgr: 

HG71 

210.6260 

-45.9354 

7.27 

1.99 

-8.00 

nil 


342 

0.52 

c 

3 

Disk 

HG72 

214.5320 

63.8186 

9.55 

1.13 

-6.69 

i8ii 


135 

0.64 

c 

14 

Unc 

HG73 

218.1200 

-5.5217 

11.22 

2.73 

-7.14 

47l?i 

40 

287 

1.38 

c 

1 

Sgr 

HG74 

219.4930 

-4.2267 

10.87 

1.32 

-7.91 

58lf 

77 

262 

0.56 

c 

17 

Sgr 

HG75 

221.6300 

-0.9168 

11.04 

1.34 

-7.61 

54lf 

45 

227 

1.43 

c 

17 

Sgr 

HG76 

221.6840 

5.2108 

11.34 

2.27 

-7.02 

47111 


232 

1.06 

c 

19 

Sgr 

HG77 

222.4390 

1.4490 

12.60 

0.90 

-5.95 

52lf 


88 

0.16 

c 

2 

Sgr 

HG78 

224.3620 

5.2677 

11.43 

2.79 

-6.53 

3911? 

4 

218 

1.34 

c 

17 

Sgr 

HG79 

225.2790 

-5.5283 

11.47 

2.06 

-7.36 

59111 

72 

257 

0.86 

c 

16 

Sgr 

HG80 

226.2300 

35.7993 

9.68 

2.31 

-7.59 

29tV 


309 

1.70 

c 

12 

HAG: 

HG81 

228.4230 

-3.8967 

11.27 

1.38 

-7.52 

5811? 

58 

222 

0.74 

c 

1 

Sgr 

HG82 

228.7960 

-13.5412 

10.76 

1.75 

-7.43 

1 + 

81 

242 

1.72 

c 

16 

Sgr 

HG83 

230.6520 

-6.4262 

11.39 

1.33 

-7.30 

rK. + 12 

oo_io 

104 

196 

0.85 

c 

1 

Sgr 
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Table 1. Continued 


ID 

RA 

Dec 

Kso 

(J-Ks)o 

MKso 

R© 

Vgsa 

Period 

Amp 

Seq 

Ref 

Feature 


degs 

degs 




kpc 

km/s 

days 

mags 




HG84 

230.6850 

-12.6305 

11.38 

1.58 

-7.47 

60t9° 

106 

234 

0.37 

C 

17 

Sgr 

HG85 

231.4140 

42.4141 

8.17 

1.95 

-7.92 

utl 

64 

322 

2.40 

C 

1 

HAG: 

HG86 

231.8480 

4.4744 

8.11 

2.01 

-7.75 

1511 

169 

303 

0.90 

C 

17 

Unc 

HG87 

239.6760 

18.8796 

10.93 

1.31 

-7.22 

4316 

164 

186 

0.35 

C 

16 

Sgr: 

HG88 

240.4270 

-12.8290 

12.32 

1.95 

-6.80 

69l^? 


193 

0.91 

C 

10 

Sgr: 

HG89 

249.1320 

-3.3937 

5.98 

3.75 

-7.35 

5li 

-40 

467 

1.89 

C 

19 

Disk 

HG90 

258.7400 

42.1733 

11.16 

1.35 

-7.67 

5911? 


233 

0.26 

C 

13 

Unc 

HG91 

268.5050 

26.4534 

9.51 

2.34 

-6.64 

1711 


199 

1.45 

C 

3 

HAG: 

HG92 

273.3720 

45.5220 

6.70 

3.80 

-6.94 

5li 

-78 

370 

2.22 

C 

19 

Disk 

HG93 

294.2910 

-35.5042 

10.14 

1.61 

-7.41 

341® 

157 

112 

0.16 

C’ 

17 

Sgr 

HG94 

294.3920 

-35.5440 

9.04 

2.02 

-8.07 

2711 

173 

367 

0.61 

C 

17 

Sgr 

HG95 

295.3690 

-32.5606 

9.39 

1.37 

-7.84 

281® 


260 

0.32 

C 

12 

Sgr 

HG96 

295.5790 

-35.3270 

9.96 

2.51 

-6.78 

231^ 

156 

231 

1.58 

C 

16 

Sgr 

HG97 

295.5890 

-32.1840 

9.94 

1.92 

-7.24 

281® 

163 

234 

1.12 

C 

16 

Sgr 

HG98 

297.2110 

-30.9750 

10.16 

3.32 

-7.01 

2811? 

171 

333 

1.62 

C 

16 

Sgr 

HG99 

298.3760 

-38.6000 

9.24 

2.00 

-7.39 

2llt 

175 

256 

1.53 

C 

16 

Sgr 

HGIOO 

300.2200 

-34.8657 

9.77 

1.32 

-7.52 

291® 


217 

1.31 

C 

3 

Sgr 

HGlOl 

300.4330 

-30.4129 

10.32 

1.25 

-7.00 

soil 


166 

0.88 

C 

20 

Sgr 

HG102 

300.7660 

-19.8179 

9.06 

1.35 

-7.15 

18li 


186 

1.48 

C 

19 

HAG 

HG103 

305.0020 

-5.5973 

8.59 

3.40 

-7.51 

1716 


431 

1.73 

C 

19 

Disk 

HG104 

305.1150 

-14.8243 

8.68 

3.11 

-7.32 

1611 

160 

356 

1.78 

c 

16 

HAG 

HG105 

308.4490 

-46.6057 

8.91 

1.57 

-7.58 

20li 


241 

0.93 

c 

20 

Sgr 

HG106 

315.0780 

-6.1150 

12.06 

3.21 

-5.39 

3111? 


145 

1.45 

c 

10 

Sgr: 

HG107 

321.5260 

-10.2295 

10.07 

0.79 

-8.41 

50l^ 


280 

0.13 

c 

14 

Unc 

HG108 

321.8190 

-30.8660 

9.48 

1.58 

-7.50 

2516 


233 

1.62 

c 

3 

Sgr: 

HG109 

328.8620 

23.7038 

5.48 

3.37 

-7.22 

4l? 


370 

1.16 

c 

18 

Disk 

HGllO 

331.3110 

0.1457 

9.28 

2.22 

-6.94 

18lt 

115 

221 

1.32 

c 

16 

HAG 

HGlll 

331.7240 

-25.1080 

8.94 

1.97 

-7.95 

241® 

70 

330 

1.60 

c 

16 

Sgr 

HG112 

333.9050 

-0.0498 

9.33 

1.48 

-7.42 

23lt 

102 

217 

0.19 

c 

1 

HAG: 

HG113 

334.2910 

-26.1176 

8.87 

2.15 

-7.68 

2llt 

64 

307 

1.50 

c 

16 

Sgr 

HG114 

335.7550 

22.2824 

9.58 

2.07 

-7.83 

3lli 


322 

1.39 

c 

18 

Unc 

HG115 

336.5810 

26.0606 

4.51 

2.71 

-7.51 

3ll 

226 

341 

1.38 

c 

1 

Disk 

HG116 

337.0440 

-13.7730 

9.88 

1.96 

-7.68 

33t8° 

-15 

288 

1.32 

c 

1 

Unc 

HG117 

340.9600 

-57.0232 

7.99 

1.54 

-7.88 

15li 

62 

274 

0.60 

c 

17 

Sgr 

HG118 

341.6190 

-27.4490 

13.27 

1.58 

-8.12 

190li 

45 

314 

0.35 

c 

17 

Unc 

HG119 

347.1460 

40.5928 

7.76 

2.60 

-7.69 

13lt 


361 

1.59 

c 

3 

Disk 

HG120 

349.8980 

-18.9401 

9.96 

1.52 

-6.60 

2lli 

37 

148 

0.89 

c 

16 

Sgr 

HG121 

350.5460 

-15.3045 

9.36 

1.49 

-6.74 

1711 


157 

0.88 

c 

11 

Sgr 


The references are numbered as follows: 1. [Totten &: Irwin] ||1998|l, 2. |Green| (|201^, 3. [Reid et al.| ||200^, 4. [Liebert et 
|al.[ <[2000l l , 5. [Lee Chen[ l [2009[|, 6. [Cruz et al.j ||2003[ |, 7. [Moody et al.[ < [1997|l ,~ Groenewegen, Oudmaijer, & Ludwig 


(|1997[|) 9. [Totten[||1998[|, 10. [Kirkpatrick et al. 


Brunzendorf|l[2001[|, 14![Gigoyan Sz, Mickaelian 


12008 , 11. 
20121, 15. 


Christlieb et al.] ( 2001[|, 12. [Cruz et al.[ ( [2007[l, 13. [Meusinger fc 


Gizis[([2002|| , 16. Mauron et al.[||2004[| , 17.[Mauron, Kendall, &: 


Gigoyanj ( [2005| l, 18. [Mauron, Gigoyan, Sz Kendall[ < [2007b| |, 19. [Mauron| l |2008[ |, 20. [Mauron et aT[ l |2014[ | 
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Table 2. Catalogue of second sample carbon stars. As for Table 1, the coordinates are given in decimal degrees (J2000), 
the (reddening corrected) NIR magnitudes and colour are from 2MASS data. Errors on distances are -^30% 


ID 

RA 

degs 

Dec 

degs 

Kso 

(J-Ks)o 

MKso 

Distance 

kpc 

^GSR 

km/s 

Reference 

HG2-1 

2.0708 

-75.6242 

11.60 

2.18 

-7.53 

67 


12 

HG2-2 

6.9679 

-70.8753 

11.04 

2.27 

-7.49 

51 


12 

HG2-3 

9.1348 

-22.9142 

14.14 

1.32 

-6.98 

167 


19 

HG2-4 

11.0646 

-71.2611 

11.36 

2.02 

-7.59 

62 


12 

HG2-5 

11.5033 

-75.3533 

11.29 

1.83 

-7.67 

62 


12 

HG2-6 

16.9675 

-69.3600 

10.97 

2.18 

-7.53 

50 


12 

HG2-7 

20.4304 

-74.0628 

11.08 

2.60 

-7.35 

48 


12 

HG2-8 

21.6450 

-70.6631 

10.34 

2.43 

-7.42 

36 


12 

HG2-9 

25.2263 

-72.1475 

10.30 

1.87 

-7.66 

39 


3 

HG2-10 

27.0325 

-71.9225 

10.31 

2.43 

-7.42 

35 


12 

HG2-11 

29.1083 

51.4225 

7.60 

1.60 

-7.40 

10 


12 

HG2-12 

33.4667 

46.6868 

8.02 

1.28 

-6.91 

10 


14 

HG2-13 

37.1960 

26.7951 

8.59 

1.47 

-7.21 

14 

50 

1 

HG2-14 

37.8151 

-2.7185 

11.52 

0.96 

-6.42 

39 

-64 

1 

HG2-15 

51.8808 

39.0806 

7.00 

1.53 

-7.30 

7 


14 

HG2-16 

56.6607 

75.7879 

8.78 

1.46 

-7.19 

16 


3 

HG2-17 

57.5445 

26.0841 

8.75 

1.37 

-7.05 

14 


18 

HG2-18 

60.4306 

8.7029 

8.11 

1.45 

-7.18 

11 


19 

HG2-19 

60.6383 

-68.2731 

11.00 

2.43 

-7.42 

48 


6 

HG2-20 

67.8925 

-72.5356 

11.16 

1.92 

-7.63 

57 


6 

HG2-21 

74.0112 

9.3719 

8.55 

1.52 

-7.29 

15 


18 

HG2-22 

74.9562 

-75.1550 

10.25 

2.30 

-7.48 

35 


6 

HG2-23 

91.6438 

73.1742 

7.84 

1.44 

-7.16 

10 


18 

HG2-24 

101.5750 

54.5260 

8.05 

1.37 

-7.06 

10 


18 

HG2-25 

105.2837 

-68.3150 

10.51 

2.22 

-7.51 

40 


6 

HG2-26 

105.9281 

23.9374 

11.67 

2.47 

-7.41 

65 


4 

HG2-27 

106.5324 

40.2004 

7.05 

1.57 

-7.36 

8 


14 

HG2-28 

118.1038 

4.5663 

7.71 

1.51 

-7.26 

10 


3 

HG2-29 

120.9702 

36.7454 

8.80 

1.20 

-6.78 

13 


14 

HG2-30 

129.3208 

-13.3539 

7.69 

1.51 

-7.26 

10 


6 

HG2-31 

134.9821 

-77.8848 

10.50 

2.23 

-7.51 

40 

69 

17 

HG2-32 

144.7517 

-10.3629 

10.71 

1.75 

-7.63 

47 

78 

1 

HG2-33 

155.5614 

-11.8608 

8.47 

1.49 

-7.24 

14 

-51 

1 

HG2-34 

175.4265 

-33.6926 

12.16 

1.64 

-7.47 

84 

-49 

17 

HG2-35 

185.7478 

21.0906 

9.97 

1.56 

-7.35 

29 

19 

7 

HG2-36 

190.9055 

2.3584 

11.44 

1.53 

-7.29 

56 

-29 

1 

HG2-37 

194.2376 

-11.7719 

9.39 

1.31 

-6.95 

19 

-44 

1 

HG2-38 

195.3270 

0.4974 

13.51 

1.62 

-7.44 

154 

-12 

19 

HG2-39 

205.6118 

-7.2565 

9.23 

1.27 

-6.89 

17 

-220 

9 

HG2-40 

212.1464 

5.1075 

7.25 

1.60 

-7.41 

9 

-37 

1 

HG2-41 

221.2036 

-1.1824 

10.74 

1.52 

-7.28 

40 

26 

1 

HG2-42 

223.4306 

-13.2182 

11.19 

2.12 

-7.55 

56 

72 

1 

HG2-43 

225.2789 

-5.5274 

11.47 

2.06 

-7.58 

65 


10 









32 A. P. Huxor and E. K. Grebel 


Table 2. continued. 


ID 

RA 

degs 

Dec 

degs 

Kso 

(J-Ks)o 

MKso 

Distance 

kpc 

km/s 

Reference 

HG2-44 

229.6677 

14.9842 

7.33 

1.44 

-7.15 

8 


20 

HG2-45 

237.9550 

-7.8469 

13.49 

3.27 

-7.07 

130 


10 

HG2-46 

245.4011 

-8.8886 

7.74 

1.39 

-7.08 

9 


3 

HG2-47 

253.1154 

-16.1274 

11.30 

1.31 

-6.96 

45 

142 

17 

HG2-48 

256.7116 

40.2096 

9.67 

1.25 

-6.86 

20 


14 

HG2-49 

259.4814 

4.6020 

10.17 

1.31 

-6.96 

27 


18 

HG2-50 

261.4760 

3.0070 

13.37 

1.33 

-7.00 

118 

32 

19 

HG2-51 

262.1073 

70.1418 

9.01 

2.50 

-7.39 

19 

37 

16 

HG2-52 

269.5650 

22.5980 

8.05 

1.41 

-7.11 

11 


18 

HG2-53 

281.7095 

-56.2341 

11.01 

1.64 

-7.47 

50 

5 

17 

HG2-54 

282.4620 

62.2910 

7.48 

1.29 

-6.92 

8 


19 

HG2-55 

285.2664 

87.0623 

7.26 

1.12 

-6.67 

6 


14 

HG2-56 

288.6007 

-78.3782 

7.37 

2.14 

-7.54 

10 

37 

17 

HG2-57 

292.9105 

-30.0417 

10.25 

1.63 

-7.46 

35 

168 

17 

HG2-58 

293.3575 

52.3578 

9.30 

2.42 

-7.43 

22 


3 

HG2-59 

293.8285 

54.6649 

8.22 

1.53 

-7.29 

13 


3 

HG2-60 

294.8768 

75.6946 

7.76 

1.61 

-7.42 

11 


18 

HG2-61 

299.6685 

77.7570 

8.65 

3.30 

-7.07 

14 


18 

HG2-62 

301.1035 

-30.1142 

10.78 

1.54 

-7.31 

41 


20 

HG2-63 

303.3270 

-23.6980 

8.08 

1.38 

-7.06 

11 

-106 

16 

HG2-64 

303.9993 

76.5857 

7.91 

1.49 

-7.23 

11 


18 

HG2-65 

304.5785 

-66.8494 

8.98 

1.80 

-7.68 

22 


20 

HG2-66 

312.0750 

10.4440 

7.81 

2.20 

-7.52 

12 


20 

HG2-67 

313.7274 

-28.4824 

10.81 

1.51 

-7.27 

41 

111 

16 

HG2-68 

329.9533 

44.3110 

7.68 

1.90 

-7.64 

12 


14 

HG2-69 

333.9954 

42.3795 

7.74 

5.35 

-6.21 

6 


14 

HG2-70 

335.4107 

33.5997 

7.22 

2.72 

-7.30 

8 


14 

HG2-71 

337.9708 

19.2946 

9.20 

1.38 

-7.06 

18 

-157 

1 

HG2-72 

341.9084 

-78.4548 

8.17 

1.75 

-7.63 

14 

237 

17 

HG2-73 

343.1504 

47.6903 

7.92 

1.56 

-7.34 

11 


12 

HG2-74 

349.3378 

-24.1952 

12.26 

1.47 

-7.20 

78 

41 

16 

HG2-75 

351.3807 

-30.1822 

10.99 

2.44 

-7.42 

48 

114 

16 
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Table 3. Variable carbon stars in Galactic satellites. Distance to the hosts are Sculptor: 86 kpc, Fornax: 147 kpc, Carina: 
105 kpc, Leo II: 233 kpc, Sextans: 86 and Sgr: 26 kpc ( |McConnachie||2012[ |. As for Table 1, the coordinates are given in 
decimal degrees (J2000), the (reddening corrected) NIR magnitudes and colour are from 2MASS data. The periods and 
amplitudes (in Catalina V band) of variability are derived from fits using Period04. 


ID 

RA 

degs 

Dec 

degs 

Kso 

(J-K)so 

MKso 

R© 

kpc 

Period 

days 

Amp 

Seq. 

Host 

SI 

14.9736 

-33.6418 

11.60 

3.24 

-8.06 

85 

524 

1.90 

C 

Scu 

S2 

14.9956 

-33.4764 

12.27 

1.13 

-7.46 

89 

195 

0.96 

C 

Scu 

S3 

39.7211 

-34.8224 

13.45 

2.02 

-7.30 

141 

242 

1.09 

C 

For 

S4 

39.7375 

-34.7760 

13.08 

1.71 

-7.78 

149 

277 

0.46 

C 

For 

S5 

39.7377 

-34.7969 

13.56 

1.26 

-7.31 

150 

189 

0.30 

C 

For 

S6 

39.8138 

-34.2523 

13.62 

1.18 

-7.44 

163 

196 

0.95 

C 

For 

S7 

40.0039 

-34.3787 

13.00 

1.52 

-7.84 

148 

128 

0.18 

C’ 

For 

S8 

40.0056 

-34.3386 

13.34 

1.35 

-7.66 

159 

231 

0.47 

C 

For 

S9 

40.0114 

-34.5303 

13.18 

1.61 

-8.00 

172 

296 

0.27 

C 

For 

SIO 

40.0278 

-34.3895 

12.62 

1.87 

-8.10 

139 

340 

0.73 

C 

For 

Sll 

40.0426 

-34.5561 

12.55 

1.52 

-8.31 

148 

334 

1.37 

C 

For 

S12 

40.0651 

-34.5675 

13.23 

1.70 

-7.81 

162 

280 

0.46 

C 

For 

S13 

40.0741 

-34.4599 

13.18 

2.24 

-7.39 

130 

273 

0.96 

C 

For 

S14 

40.0831 

-34.5527 

13.91 

1.36 

-7.24 

170 

190 

0.83 

C 

For 

S15 

40.1043 

-34.4829 

13.51 

1.32 

-7.64 

169 

226 

0.45 

C 

For 

S16 

40.1301 

-34.4789 

13.07 

1.67 

-7.79 

149 

275 

0.76 

C 

For 

S17 

40.2177 

-34.6231 

13.66 

1.73 

-7.29 

155 

219 

0.36 

C 

For 

S18 

40.2223 

-34.2037 

13.26 

1.77 

-7.44 

138 

239 

0.85 

C 

For 

S19 

40.2648 

-34.8015 

12.69 

1.75 

-8.03 

139 

315 

1.19 

C 

For 

S20 

40.2847 

-34.5978 

13.43 

1.59 

-7.41 

147 

222 

0.31 

c 

For 

S21 

100.3062 

-50.9068 

12.66 

1.27 

-7.21 

94 

181 

0.17 

c 

Car 

S22 

100.4230 

-50.9689 

12.34 

1.40 

-8.05 

119 

284 

0.66 

c 

Car 

S23 

100.4353 

-51.0058 

12.69 

1.55 

-7.60 

115 

241 

0.41 

c 

Car 

S24 

100.5431 

-50.9400 

12.78 

1.15 

-7.42 

110 

193 

0.38 

c 

Car 

S25 

153.8581 

-2.0757 

11.95 

2.10 

-7.99 

97 

348 

1.41 

c 

Sex 

S26 

168.3034 

22.1872 

14.79 

1.44 

-6.89 

217 

164 

0.42 

c 

LeoII 

S27 

168.3359 

22.1878 

14.81 

1.56 

-7.03 

234 

183 

0.82 

c 

LeoII 

S28 

288.9922 

-30.5467 

9.73 

1.94 

-7.66 

30 

281 

0.47 

c 

Sgr 

S29 

289.7724 

-30.9072 

9.16 

1.63 

-7.54 

22 

240 

1.88 

c 

Sgr 

S30 

290.9395 

-34.5209 

9.74 

1.28 

-7.89 

34 

252 

1.26 

c 

Sgr 

S31 

292.2203 

-29.5899 

10.16 

1.42 

-6.85 

25 

160 

0.55 

c 

Sgr 

S32 

295.7120 

-33.5759 

10.14 

1.33 

-7.76 

38 

241 

0.92 

c 

Sgr 







